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ABSTRACT 

The design of a class of special purpose computing machines 
which compute by counting is systematically developed. The basis 
of the design philosophy is to limit the basic building elements to 
three fundamental units and to develop the method of synthesis such 
that these three building elements are represented as operational 
units. In particular, the three basic building elements are (l) the 
binary rate multiplier which is a means of scaling down a pulse 
stream to some specified fraction, (2) the counter, and (3) the 
anti -coincidence circuit which is a means of separating pulses 
arriving at the counter simultaneously. The computational errors; 
i.e., rounding -off error and truncation error, introduced into the 
machines when these elements are treated as operational ’units are 
studied in detail. The method of synthesis is explicitly stated 
and a wide variety of machines obtained directly from this synthe- 
sis are presented. Finally, a series of machines is presented for 
interpolation and extrapolation of a function which is available 
only as empirical data. 
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CHAPTER I 


PRINCIPLE DESIGN ELEMENTS 
Introduction 

Computers are usually divided into two broad categories, analog 
and digital. Analog computers represent variables as physical quan- 
tities . The solution of a problem in an analog computer is attained 
by constraining a physical model of the problem to be solved, and 
measuring the variables . The ability to program a wide variety of 
problems is achieved by having functional components available (e.g., 
adders, multipliers, integrators) and interconnecting them by means 
of a patchboard. The resulting interconnection is scaled to matcn 
the desired equation. On the other hand, digital computers represent 
variables as discrete quantities. The usual method of solution of a 

J 

problem in a digital computer is attained by sequencing a sequenc • of 
instructions through the fetch-execute cycle of its control unit. 
Another class . f digital computers, known as incremental computers, 
combine the parallel functional simplicity arid speed of analog com- 
puters with the ability of attaining computational precision which is 
not dependent on precision of measurements. Such computers attain a 
speed advantage over conventional general purpose ccjpuiers by trans- 
mitting and processing only partial words in a number of parallel 
arithmetic organs rather than the whole words needed by the fetch-, 
execute cycle. Moreover the digital nature of these computers permit 
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the problem solution to he repeated exactly and therefore does not 
possess the drift characteristic of analog computers . Beyond a doubt 
the incremental computer which has found the most interest in the 
literature is the digital differential analyzer; i.e., DM. This 
computer can be viewed as a digital analogy of an analog computer. 

The usual design practice in each of these machines is to permit them 
to solve a large spectrum of problems. When a computer need arises 
for a special purpose application, this versatility is felt as a cost 
factor. 

A class of incremental techniques which has been used in real 
time control is a class known in the industry as countup-countdown 
techniques . The basis of these techniques is to represent data by a 
unitary code. For example, the number 28 is represented by 28 pulses. 
A function may be represented by counting the sequence of pulses in a 
foruard-b&ckvard counter or converting them directly into an analog 
quantity (e.g., by a stepping motor) for analog processing. Conse- 
quently, when a real time application deals with continuous-smoothly 
varying functions, countup-countdown techniques offer a simplicity 
and economy of hardware which is hard to beat with computing systems 
designed to handle a large spectrum of problems. 

The purpose of this thesis is to investigate countup-countdown 
techniques with the objective of demonstrating that they can, in 

r 

fact, be used to generate a wide variety of non-trival functions. 

This will be done by displaying a circuit which will generate each 
function. However, since the techniques upon which we base this 
thesis are described in the literature only in an ad-hoc manner 
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(Refs. 4, 5, 9, 10, 12, and 19), we wi'J. be specific as to which cir- 
cuits we will permit as basic building elements. In particular, the 
fundamental units which we will permit are (l) the binary rate multi- 
plier (abbreviated BFM) which is a means of scaling down a pulse 
stream to some specified fraction, (2) the counter, and (3) the anti- 
coincidence circuit which is a means of separating pulses arriving 
at the counter simultaneously. In order to strengthen our argument 
we will avoid completely the explicit use of adders and subtractors. 
A succinct recapitulation of the purpose of this study is to system- 
atically develop and demonstrate the versatility of techniques based 
on counting for solving sophisticated and practical special purpose 
computer design problems. 

Our method of synthesis will be to describe the principle 
building elements as operational units and then proceed by opera- 
tional techniques to show how to fabricate the various machines. In 
particular, a first order difference equation can be represented by 
a counter, and approximate integration can be attained by using a 
counter in cascade with a binary rate multiplier. These principle 
design elements are described in this chapter. 

It is to be expected that the results obtained by operational 
means will deviate from the actual results due to the finiteness of 
the machine and the approximation implied by our synthesis. A dis- 
cussion of these approximations is presented in CHAPTER II. This 
chapter is supplemented by Appendices A and B where some quantitative 
results are presented related to the computational accuracy of the 
BRM. In CHAPTER III we explicitly state the method of synthesis and 
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demonstrate it by deriving a wide variety of representative machines. 
Some of these machines have been simulated on a general purpose com- 
puter and these results are also presented and discussed in CHAP- 
TER III. In CHAPTER IV the specific problems of constructing poly- 
nomial generating machines are considered. In particular. * family 
of machines are given for interpolating and extrapolating values of a 
function defined only by empirical data. 

Binary Rate Multiplier 

A binary rate multiplier (abbreviated ERM) is a means of scaling 
down a pulse stream to some specified fraction. A logic diagram of 
a BRM which is built out of standard logic elements is shown in 
Fig. 1.1a. This circuit is described in detail in several of the 
references (e.g., Refs. 4 and 10). Consequently, a brief description 
will serve our purposes. The input pulse stream is applied directly 
to the binary counter whose value is denoted by . . . X£X^. 

Each flip-flop of the counter is operated as a trigger. For every 
two input pulses to a trigger two output Folses are produced; one 
pulse when the flip-flop makes a 0 to 1 transition called an a pulse 
and one when the flip-flop makes a 1 to 0 transition called a 
R pulse. The p pulse is used to trigger the successive stage of the 
counter. The a pulses are gated through gated pulse generators and 
mixed through a NOR element to produce the desired fraction of the 
input pulses. This simple mixing technique may be used because the 
a pulses from the various stages are separated in time from each 
other. This timing factor is shown in Fig. 1.1b. 

The quantitative relationship of a ERM may be expressed as 




<b) Timing diagram. 

Figure 1. 1. •• Binary rate multiplier. 
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follows: If Ax is the number cf input pulses, the number of out- 

put pulses produced by the k tb stage of the counter is Ax • 2"^, 
This multiplicative relation will remain valid over any interval for 

T/- 

which Ax is a multiple of 2 pulses. If y_ k is the level set- 
ting of che k stage gated pulse generator, the number of output 
pulses which may be gated through this stage will be y_^ Ax • 2 . 

Since the output pulses from the various stages are simply mixed, 
the number of output pulses Az of an n stage BRM over any inter- 
val Ax which is a multiple of 2 n pulses will be the sum of all 

the nulses gated through all the stages. This output is 

n 

Az = Ax y wi 2 -i (l.l) 

i=l 

The quantity y = t y-i 2 ' 1 is a binar y number. Therefore 

fcl 


Eq,. (l.l) may be written as 


Az = y Ax 

(1.2) 

where the range of y is 


0 < y < 1 - 2”* 1 in steps of 2 r ' 

(1.3) 


If y remains constant over a Ax interval of 2 n pulses, 
then the output shown by Eq. (1.2) remains exact. However if Ax. is 
less than 2 n pulses then this multiplicative relationship remains 
valid only on the average. This can be demonstrated as follows: If 

Ax is the number of input pulses into e. n-stage BRM starting with 
counter value x, and whose gated pulse generators are set to value 
y, then the output for this machine is Az x . Since there are 2 n 
possible starting values, t.hen there are 2 n possible different 



machines. The average output; denote it by As, over all of these 
2 n different machines is 


n -» 


2 n As = 


Az 


x 


x=0 


(1.4) 


The total pulse output ovar all of these machines is 



Ass, 


pulses. This is equivalent to putting 2 n Ax successive input 
pulses into a single machine since each transition over all of these 
machines is attained Ax times. For example, the transition ending 
with counter value x is attained by the Ax machines starting out 
with the counter value prior to x. Therefore, the total number of 
output pulses over all of these machines is also given by Eq. (1.2). 

2 n -l 


Ae x = y * 2 Ax 


( 1 . 5 ) 


Combining Eqs. (1.4) and (1.5) we have 

Az = y Ax (1.6) 

Because of the approximate nature of Eq.. (1.2) when Ax is 
less than 2 n pulses, we will calculate the specific output sequence 
in demonstrating specific machines. For these calculations, the 
pulse stream shown in Fig. 1.1b may be displayed in vector form. 

This will be called the p-sequence. Each position of the vector in 
this sequence represents the possible output at a particular pulse 
time from a stage of the BRM. The p-sequence for a two, three, an<_ 
four stage BRM is displayed in Table 1.1 
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TABLE 1.1 - EXAMPLES OP p-SEQTJENCES 


Pulse 

■ 2 Stage 
p-sequence 

3 Stage 
p-sequence 

"" — 1 

4 Stage 
p-sequence 

1 

10 

100 

1000 

2 

01 

010 

0100 

3 

10 

100 

1000 

4 

00 

001 

0010 

5 


100 

1000 

6 


010 

0100 

7 


100 

1000 

8 


000 

0001 

9 



1000 

10 



0100 

11 



1000 

12 



0010 

13 



1000 

14 



0100 

15 



1000 

16 



0000 


The p-sequences given above assumes that the BRM counter start- 
ing value is zero. If another starting value is used then its 
associated p- sequence can be easily obtained. Moreover, if an in- 
terval greater than 2 n pulses is used, then the p- sequence can be 
obtained by repeating the p-sequence given above. 

Hie sequence of output pulses may be calculated by multiplying 
oit-by-bit the p-sequence with the respective values ..of the level • 
settings of the gated pulse generators . This process is illustrated 
below by two examples. 
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Example A: 


Example B: 



The first matrix, in each of these examples, is the p-sequence. 


The 


next matrix represents successive values of the gate settings. When 


these two matrices are multiplied, the result is developed along Lie 
diagonal of the resultant matrix. This result is shown as a vector 
on the right hand side. 

The expected output value of Example A by Eq. (1.6) is 35/8 
pulses for the 8 input pulses. -However, as shown by actual computa- 
tion, the ERM yields zero output pulses . On the other band, the 
expected value of Example B by Eq. (1.6) is 2l/8 pulses for the 
7 input pulses. The above computation yields 7 output pulses. Both 
of these examples are pathological cases in the use of the BKM. The 
approximate nature of Eq. (1.6) can ordinarily be expected to yield 
more realistic results. Some of these results are presented in 
CHAPTER III. 

The method of synthesis to be presented necessitates that the 
ERM operate on signed quantities . In particular, the level -setting of 
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Figure 1.2. - Signed four-stage bidirectional counter. 


the gated pulse generators and the counter input pulses must he 
signed quantities, and the BRM is to yield a signed pulse output. 

If the output pulses are accumulated in a counter then the sign of 
the pulse will determine the direction of counting. If the output 
pulses are used to drive a stepping motor, then the sign of the pulse 
will determine the direction the stepping motor is to turn. Through- 
out this discussion we will consider that the signs- --of. various quan- 
tities are available through level logic . Consequently, the output 
sign can be obtained from the input signs by an exclusive OR circuit. 

Counter 

The purpose of the counter in the machines which will be con- 
sidered are twofold; (l) to accumulate the pulses arriving at the 
counter in order to display the total. number of counts, and (2) to 
set the levels of the ERM's. In the first application the counting 
sequence can be any desired sequence for a terminal device. In many 
real time applications the output pulses may not be accumulated di- 
rectly but are converted to an analog quantity for analog process- 
ing (e.g., by a stepping motor). In the second application, the 
counting sequence must be compatible with the BRM. This general re- 
quirement can be met by the circuit displayed in Pig. 1.2. 

A number is represented in vhis counter by magnitude plus sign. 
As had been stated earlier the signs are represented by level logic. 
The counter counts down :.n magnitude when the input pulse and counter 
are opposite in sign .and counts up in magnitude when t: -e -counter and 
input pulse have the same sign. The circuit is designed sb that the- 
co-pulses are used to count down and the ft pulses used to count up. 
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There are two representations of zero; that is, minus zero and 
positive zero. When the counter value is at +00 ... 0 and a -1 
pulse arrives then the counter is set to -00 . . .01. This end cor- 
rection is accomplished in three steps. The normal sequence first 
changes the counter value to +11 ... 1. The magnitude is then 
corrected in the second step to +0 . . .01. Finally, the sign is 
changed to -00 . . .01. The sign is changed last so that the £ 
pulses generated when the magnitude is corrected do not propagate to 
the successive stages of the counter. In a similar manner to that 
just given, the counter is set to +00 . . . OX when the counter values 
is -00 ... 0 and a +1 pulse arrives. The down counting sequence 
fcr a three stage counter is given in Table 1.2. 


TABLE 1.2 - DOWN COUNTING SEQUENCE 


-1 Input pulse 

+1 Input pulse 

+111 

+13.0 

+101 

+100 

+011 

+010 

+001 

+000 

+111 -> +001 - -001 

-111 

-110 

-101 

-100 

-Oil 

-olo 

-001 

-000 

-111 -» -001 +001 


The up counting sequence utilizing the |3 pulses of the flip- 
flops. is given in Table 1.3. 
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TABLE 1.3 - UP COUNTING SEQUENCE 



+1 Input pulse j 

-000 

+000 

-001 

+001 

-010 

+010 

-Oil 

+on 

-100 

+100 

-101 

+101 

-no 

+110 

-in 

+ni 


Since the signs of both the pulse output of the BRM and counter 
value are to be processed by level logic, then the activation of the 
up-down line is accomplished by an exclusive OR circuit. This is 
obvious from Table 1.4. 

TABLE 1.4 - COUNTER SIGN CONTROL 


Sign 

input 

pulse 

Sign 

counter 

Line 

activated 

+ 

+ 

Up 


- 

Down 

- 

+ 

Down 

- 

- 

Up 


Anti -Coincidence Circuit 

Pulses arriving at a counter simultaneously must first be 
separated before they are entered into the counter. Hie circuit 
that accomplishes this task is called an anti-coincidence circuit. 
Fundamentally, this circuit necessitates storing each pulse as it 
arrives. Each stored pulse is 'then presented to the counter accord- 
ing to a fixed program. The circuit configuration which can accom- 
plish this task for two inputs is shown in Fig. 1.3. 

The operation of the circuit given in Fig. l._3 is as follows.: 
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If a pulse from input 1 exists, it i.s stored in flip-flop 1. If a 
pulse from input 2 exists, it is stored in flip-flop 2. It will be 
noted that these two inputs can arrive simultaneously. Two pulses 
are emitted by the clock which are separated from each other. If 
flip-flop 1 had been set by input 1, it is reset by the clock 
pulse C]_, which in turn generates an output pulse. If flip-flop 1 
had not been set, no output pulse will appear in the output. Flip- 
flop 2 is reset and an output is similar 1> generated by clock 
pulse Cg. Since clock pulses and Cg are separated, the 

corresponding output pulses are also separated. 

Since the sign of a pulse is processed by level logic, the sign 
need not be stored before they are presented to the anti -coincidence 
circuit. However, when a pulse is presented to the counter, its 
sign must also be presented. This may be simply accomplished by 
shifting the sign level to a flip-flop by the separated clock pulses 
Ci and C2. This circuit is also shown in Fig. 1.3 where S’s and 
S's are the sign levels of the pulses and their complements, respec- 
tively. 

If more than two inputs arrive at the counter simultaneously, 
then a need arises for a circuit other than a simple clock to sepa- 
rate the stored pulses A simple binary counting sequence such as 
the leftmost sequence shown in Table 1.5 will serve this purpose. 
However, it will be noted that while this sequence can generate more 
than two steps, the a and 0 pulses from the various flip-flops are 
net separated, so consequently can not both be used. 

All the sides of the flip-flops could be used if the counting 
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sequence utilizes ~ unit distance code. Such a code would guarantee 
that not more than one flip-flop would change state at any step of 
the counting sequence. Consider the Gray code counting sequence 
given by the middle sequence in Table 1.5. This counting sequence can 
be used to separate as many as six inputs arriving at the counter 
simultaneously. However, this requires the counter itself to go 
through eight steps. An example of a counting sequence which may be 
used to handle six inpats and yet go through only six steps in the 
counting sequence is given by the rightmost sequence in Table 1.5. 


TABLE 1.5 - PULSES GENERATED BY SEVERAL COURTING SEQUENCES 


Counting 

sequence 

Pulses 

generated 

Counting 

sequence 

Pulses 

generated 



Counting 

sequence 

Pulses 

generated 

000 

— a 

000 

— a 

000 

— a 

001 


001 

-a- 

001 

-a- 

010 

—a 

on 

--P 

on 

a — 

on 


010 

a — 

111 

-0 

100 

--a 

no 

— a 

no 

-0- 

101 

-a0 

in 

-0“ 

100 

0— 

no 

— a 

101 

"0 



111 

000 

100 

0" 




Schematic Representation 

The three circuits described In this chapter are the principle 
design elements. However, in describing the machines promised by 
this thesis, these three circuits will be represented as operational 
units. The advantages to be gained by using operational units rather 
than these circuits are twofold. First, the method of synthesis can 
be more clearly presented. Secondly, a considerable hardware reduc- 
tion can usually be realized when the composite machine is con- 
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sidered. These simplifications arise when all of the features of 
these basic circuits are not required. A checklist of the features 
which when removed would simplify the basic circuits would include: 
(l) sign control of BRM may not be required, (2) counter may not 
be required to both countup and countdown, (3) two ERM's may receive 
the same input pulse stream with the result that one BRM counter may 
be used with two sets of gated pulse generators, and (4) the level 
setting of the BRM may be constant with the result that a scaling 
circuit (see Ref. 10) rather than a BRM may be used. These three 
circuits are, however, sufficient as principle design elements. 

The three principle design elements as operational units arc 
presented in Pig. 1.4. The BRM is represented by the schematic dia- 
gram shown in Fig. 1.4a. The value y in this diagram is less than 
one and is obtained from level logic. The quantities Ax and ts. 
are the input and output pulse streams, respectively. The input- 
output relation for this diagram is expressed 'ey Eq. (1.2). Alter- 
nately, the BRM is represented by the schematic diagram shown in 
Fig. 1.4b when the value of y remains constant, in these cases 
the BRM may be replaced by a scaling circuit in the final design. 

Fig. 1.4c presents the schematic diagram of the counter. The quan- 
tity is an input pulse stream and z is the output which may be 

used in level logic. When the counter is used to set the levels of 
the gated pulse generators of a BRM a scale reduction of 2“ n is 
implied by the connection. At times this scale reduction will be 
shown explicitly by tt same diagram shown in Fig. 1.4b. The initial 
conditions of a counter may be shown explicitly by inserting it in 
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the box; i.e. , z (o)‘ Considering the counter value as a function of 
iterative steps then the input-output relation may be expressed by 
the first order difference equation 

z (k) = z (k-l) + ^(k-l) 

The value of z/^) in Eq. (1.7) in terms of the initial condition 
of the counter is 


z (k) “ z (o) 



( 1 . 8 ) 


Fig. 1.4d represents the schematic diagram of an anti -coincidence 
circuit. This circuit accepts multiple pulse inputs and produces 
a single pulse output. The design of this circuit is such as to 
permit the input pulses to arrive simultaneously. However, at times 
it will be convenient to use this schematic diagram for multiple 
pulse inputs even if the pulses are known to be separated. This 
usage, therefore, should permit a corresponding simplification in 
the final design. 



IS 
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CHAPTER II 

ANALYSIS OF MACHINE COMPUTATIONAL ERRORS 
Classification of Errors 

The difference between the actual outpux of a system and that 
given by a theoretical model is considered the error of the system. 
If f represents the theoretical process given by the model, x 
represents the values of the arguments, f a represents the actual 
process, and x a represents the values of the arguments vitiated 
by previous calculations, then the total error £ is given by 

£ = f(x) - f a (xa) (2.1) 

It is convenient to subdivide the error into the error propa- 
gated from previous calculations and error generated locally. The 
sum of these two errors is also equal to the total error as is evi- 
denced by rewriting Eq. (2.1) as: 

£ = f(x) - f(x a ) + f(x a ) - f a (^) (2.2) 

The quantity f(x) - f (:% ) is error propagated from previous calcu- 
lations and is called the propagated error. The difference between 
the value calculated locally b^ the model and the value generated by 
the actual process; i.e., f(:%) - f a (xa,), is called the generated 
error . 

von Neuman and Goldstine (Ref. l) classified the generated 
errors into four categories according to their source. These four 



sources are simply listed as follows: 

(1) Model errors 

(2) Input errors 

(3) Truncation errors 

(4) Rounding-off errors 

The last two errors are of primary concern to the numerical 
analyst and are called computational errors. Truncation errors re- 
sult from expressing transcendental operations as numerical pro- 
cesses. For example, if a transcendental function is evaluated by 
an infinite series or as the fixed point of a process, then the 
truncation error is the error introduced by terminating the evalua- 
tion short of the limit goals.. In the case of an iterative process, 
this error is called the iterative truncation error. Rounding-off 
error is introduced into the resultant after each arithmetic opera- 
tion. In conventional digital computers, this error is introduced by 
rounding or chopping a number such that it can be represented by a 
register of fixed length. It can be viewed as an error in the 
arithmetic processes. 

If f c represents the numerical approximation of the theoreti- 
cal process f, then the difference fC^) - f c (x^) is the generated 
truncation error and f c (xa) - f a (x a ) is the generated rounding 
error. The generated error is equal to the sum of these two errors 
as is evidenced by 

* f c (x a^ ’ f a (x a^ 
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Computational Errors of .<RM 

Fig. 2.1 represents a counter in cascade with a BRM. If y^) 
in this figure represents successive values of y, each of which re- 
main constant over some interval of 2 n pulses, then the value of 
the counter can be expressed by the following difference equation. 

z (k) = z (k-l) + y(k-l) ** ( 2 - 4 ) 

If z(q) represents the initial value of the z counter, then 
Eq. (2.4) may be expressed as 

V •» 

z (k) = z (o) + IT y(i) < 2 * 5 ) 

1=0 

This equation is recognized as Euler's (rectangular) integration. 

A model of this process which is convenient for machine synthe- . 
sis is presented in Fig. 2.2. The deviation of the results given -by 
the model from that given by Eq. (2.4) is the truncation error. 

The counter z in Fig. 2.1 may be viewed as a lower register 
z 1 consisting of n stages and an upper register z u of an arbi- 
trary number of stages. Thus, after the first iteration of 
Eq. (2.4), z^ contains the fraction y(Qi ? n of input pulses. 
Rounding of the upper register may be accomplished by presetting 
z^ to one -half of the maximum counts possible in z^, and chopping 
may be accomplished by presetting z-> b to zero. After the second 
iteration, y(q) 2*- 1 ' pulses are added to the counter z. As a. : result 
of this iteration, Zj may or may not overflow into z u . Proceeding 
^ih' this manner, it is noted that z u represents- the single -precision 
rounded or chopped sum shown: in Eq. (2.5). 

The above description has been presented only in order to put 



Figure 2.2. - Integration model. 
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the rounding error in the framework of conventional computers. A 
more realistic operational procedure may be realized is Ax is a 
subinterval of less than 2 n pulses. The advantages to be gained in 
this situation are twofold; (l) the speed of the computation may be 
increased, and (2) a hardware savings may be realized by reducing or 
eliminating the lower register. The output in this case is v imitated 
also by an error in multiplication. In particular, if Ax is a 
subinterval of one pulse then the output is in error only by the 
error in multiplication. Because of the importance of Eq. (2.4) in 
this study, the multiplication error of this equation is presented 
in Appendices A and B and in the following sections. 

Multiplication Error Formulas 

Starting the BRM counter with zero the ^actual output of the ERM 
is given by . 

z = Entier (y Ax + l/2) (2.6) 

This function together with a plot of Eq. (1.2) is given in 
Fig. ; (2.3) for a three stage BRM for the various values of y. The 
difference E between these two quantities; i.e., 

E = Entier (y Ax + 1/2) - y Ax (2.7) 

is plotted for this three stage BRM in Fig. (2.4). The difference 
E, when only one stage of an n stage BRM is gated by y, can also 
be expressed systematical]^ in tabular form as shown in Table 2.1. 
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TABLE 2.1 - MULTIPLICATION ERROR E OF BRM 



An inspection of these tables shows that the error associated 


with the various stages of a BRM may be expressed more concisely in 
algebraic form as shown in Table 2.2. 

TABLE 2.2 - MULTIPLICATION ERROR E 


IN ALGEBRAIC FORM 



For an arbitrary value of y, the value of E is the linear 
combination of the values shown in Table 2.2. This bilinear form 
is shown in Eq. (2.8) for an n stage BRM. The element subscripts 
of the Boolean vectors x and y (i.e,, vectors whose elements are 
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0 or l) which are shown in this equation correspond to the stage 
numbers of the BRM shown in Fig. 1.1. 



( 2 . 8 ) 


(2.9) 


\ 




_1 

~4 

1 
. 2 


0 


1 

e 

l 

4 

1 

:2 


In the formulation of E, the maximum /a lues of the output 
of the BRM were reflected at the points of discontinuities . It 
will be observed that just prior to these points the error is one 
quanta less than that shown by E. A formulation of F in which 
the minimum values are reflected at the points of discontinuities 
can be obtained in a manner similar to that for obtaining E. The 
quantity F, when only one stage of an n stage BRM is gated by 
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y, is shown in tabular form in Table 2.3 (only a few cases are 
exhibited). 


TABLE 2.3 - MULTIPLICATION ERROR F WHEN 
ONLY ONE STAGE IS 0 .TED 



It will be observed that the F values equal the E values 
except at the points where the discontinuity occurs. At these 
points F equals -l/2 while the corresponding E value equals 
+l/2. The C values shown above correspond to the 2’s comple- 
ment of the x values. It will be observed that the F values are 
identical in terms of C to the negative values of E. Conse- 
quently, it can be asserted that F in terms of C is just the 
negative of E . 



F = -(C^Cg, . . C n )M * j 

V-n/ 


( 2 . 10 ) 
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An example will help clarify these formulas. In this example 
the value of y is 101, and the values of E and F are calcu- 


lated for successive ERM counter values. 



This example is shown in graphical form in Fig. 2.5. 

Consider the difference £ - F. For a three stage BRM this is: 


/ y -A 

E - F = (x 1 ,x 2 ,x 3 )Mj y 2 J 

V-s/ 


+ (c 1 ,c 2 ,c 3 )m| 



( 2 . 11 ) 


The premultiplier to the vector y; i.e. , 

(x 1 ,x 2 ,x 3 )M -J (C 1 ,C 2 ,C 3 )M » (x x + C 1 ,x 2 + C 2 ,x 3 + C 3 )M 
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When a BRM counter starts out with an arbitrary value , then 
the starting value must enter into the error formula as a param- 
eter. Theoe error formulas are given explicitly by Eqs. (2.12) 
and (2.13). In particular, Eqs. (2.12) and (2.13) reflect the maxi- 
mum and minimum values of the actual output at the points of dis- 
continuities, respectively. In these formulations x and Xg rep- 
resents the value and the initial value of the counter, and C 
and Cg represents the 2's complements of these values. The sub- 
scripts on those literals represent, as before, the stage of the 
BRM. In Eq. (2.13), Xgg identifies the right most counter bit 
whose value is 1; e.g., for the counter val - 100 then Xg^y_^- y_g, 
for counter value Oil then xggy_R = y_^, etc. . 


G = (xi - x sl ,x 2 - x S2 , 


H “ ~ x SR y -R “ (°1 " C S1' C 2 _ C S2> 



( 2 . 12 ) 


(2.13) 


Multiplication Error Bounds 

The maximum positive error and the minimum negative error for a 
n stage BRM whose counter starts out with zero may be obtained by 
an analysis of Eqs. (2.$) and (2.10), respectively. These values 
will then form a bound of the deviation .of the BRM from that of 
exact multiplication. This analysis is presented in Appendix A. It 
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is shown in that analysis that for an n stage ERM these values 
are: 

Emax(n) * ~ + | + (2.14) 


F„ n -„(n) = ~ - 

min' ' 1Q c 


7 n (-l) n 


18 


9 • 2 


on 


(2.3.5) 


Eq. (2.14) is plotted together with Eq. (2.15) in Pig. 2.6. As 
a hy-produet of developing Eqs. (2.14) and (2.15), it was necessary 
to find the points where these values occurred. These points are 
tabulated for various ERM in Tables 2.4 and 2.5. 

Appendix B presents an analysis of a BRM whose counter starts 
out with an arbitrary value. The basis of this analysis is to use 
Eq. (2.12) to obtai' he maximum positive error and to use 
Eq. (2.13) to obtain the minimum negative error. For a n stage 
ERM these values are: 

->n-l 


G ( n ) = ± + E - (- 1 ) 

naax' 11 ' - 


1 

9 3 


H . (n) = J. . £ . -LlL ■ 

•‘hnin' 11 ' 9 3 g . g n-l 


(2.16) 

(2.17) 


These values form a bound for the generated round-off error. 
Fortunately, these values are + a!en on at only two points of the 
BRM (for n > 2), and therefore one can expect better results than 
would be predi-: ■Jd by these values. ..These values are plotted in 
Fig. 2.7 and are also presented together with the points at which 
they occur in Tables 2.6 and 2.7. 

The problem which has been considered in Appendices A and B . 
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TABLE 2.4 - x AND y VALUES FOR E max 


n 



x 7 • • * X 1 

y-i--*y-7 

®max 

2 

11 

11 

■ 

n 

3/4 

3 

101 

101 

9 

111 

7/8 

4 

1011 

1101 


1011 

17/16 

5 

10101 

10101 

non 

11011 

39/32 

6 

101011 

110101 

110101 

101011 

89/64 

7 

1010101 

1010101 

1101011 

1101011 

199/128 | 


TABLE 2.5 - x AND y VALUES FOR F mi 


n 

BMS 


1023 

BBEBI 

F • 
A mxn 

2 

01 

11 

01 

n • 

-3/4 

3 

Oil 

101 

001 

1X1 

- 7/8 

4 

0101 

1101 

oon 

ion 

-17/16 

5 

01011 

10101 

00101 

non 

-39/32 

6 

010101 

110101 

ooion 

101011 

-89/64 

7 

0101011 

1010101 

0010101 

nonn 

-199/128 


TABLE 2 

.6 - X y Xg y 

AND y 

VALUES FOR 

^max 


n 

022 

X7. . .x x 

y-i---y -7 

X57. . *Xgl 

r — , 

Xy • * 

222 

%>ax 

2 

01 

10 

01 

00 

n 

n 

3/4 

3 

.001 

no 

on 

010 

101 

101 

9/8 

4 

0101 

1010 

0101 - 

0010 

noi 

ion 

23/16 

5 

00101 

11010 

oion 

01010 

10101 

10101 

57/32 

6 

010101 

101010 

010101 

001010 

noioi 

101011 

135/64 

B 

0010101 

noioio 

oioion 

0101010 

1010101 

1010101 

313/128 


- TABLE 2:7 - X) xg, AND y VALUES FOR H min 


n 

Xg?. . .Xgi 

X?. • .X X 

222 

x S7 . . .x sl 

mm 

■m 

ysflsyQ 

2 


01 

n 

n 

01 

11- 

3 


on 

101 

. m 

001 

in 

'4 


oon 

ion 

ion 

0101 

noi 

5 


oion 

10101 

- non 

00101 

non 

6 

noioi 

ooion 

101011 

101011 

010101 

110101 

7 

1010101 

oioion 

1010101 

1101011 

0010101 

noion 
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Figure 46. - Multiplication error 
bounds (zero starting). 
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n (number of stages) 

Figure 2. 7, - Multiplication error 
bounds (arbitrary starting). 
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and in this section is illustrated in Figs. 2.8a and 2.8b. The 
actual and desired outputs for a 2 and 3 stage BRM are plotted in 
these figures for all starting values. As is illustrated, finding 
the multiplication error bounds by graphical means is not trivial. 
The points labeled E ma y , F m i n , Graa.-, and in these simple 

cases agree with those predicted in Appendices A and B. 

Error Bounds 

The usual formulation of tne error problem is to calculate 
bounds of the total error based on bounds of the generated errors. 

In the case of Euler's integration these results are available in 
the literature (e.g. , Refs. 15 art 16). Since the analysis by which 
the bounds are obtained is based on worst case conditions, the re- 
sults are usually too pessimistic for design purpose. In particu- 
lar, P. Henrici (F.ef. 15) has presented the complete analysis of 
the initial value problem 

y' = f(x,y), y(a) = n (2. IB) 

approximated by Euler's integration. The bounds that he presents 
are in terms of Lips chi tz function; i.e., 

El ( x) = L > 0 

x L = 0 (2.19) 

where L is a constant such that for any x in the interval 
a < x < b and any "two values y ar:d y 

|f(x,y) - f(x,y*)| < L|y - y*| (2. 20) 

The bounds for the truncation error t(j c ) and the accumulated 
rounding-off error x ( k ) are given by 
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0 2 0 I 3 I 2 0 2 3 1 


x, BRM counter 
(a) Two-stage. 

Figure 2. 8. - Output of BRM for all starting values. 


icr* 












024 6013571 246 0 235713 

x, BRM counter 

(b) Three-stage. 

Figure 2. 8. - Continued. Output of BR.V for all starting values. 








46024 571 35 

x, BR<\rtc( 


(b) Concluded. 
Figure 2. 8. - Concluded. Output 
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t(k) < ^ K(x( k ))E L (x( k ) - a) (2.?-l) 

~(k) <^ E L (x (k) - a ) < 2 ‘ 22 ) 

where N(x) = 1/2 maxjy"(t)| for all t in the interval a < t < x 
and a is the maxima” 1 , local rounding error. 

The total error £/• ' is bounded by the sum of these two 
bounds; i.e., 

e (k) < ^ac 5f(x k )E L (x( k ) - a) + ^E L (x( k ) - a) (2.23) 

Consider now the generation of this function by use of a BRM. 
The function f(x,y) is used to set the level of the BRM, and the 
output of the BRM is summed in a counter which represents the value 
of y (see Eq. (2.18)). If the value of f(x,y) is updated ? very 
Ax pulses, then the bound given by Eq. (2.23) may be applied 
directly to this process. Suppose the interval size is chosen such 
that the error bound given by Eq. (2.23) is a minimum . If f(x,y) 
is updated every pulse . instead of every Ax pulses, we could, expect 
that the actual error would be smaller than this minimum. In any 
event we will take this minimum as the error bound for the function 
generated by the BRM. The minimum of the right-hand side of 
Eq. (2.23) is 

5( k ) £ VN(x(k)) E I,( x (k) " a ) (2.24) 

The value of e in Eq. (2.24) may be obtained from Eq. (2.17). 
This equation may be used to form a bound of the multiplication 
error. Moreover for large values of n, this may be approximated by 
7/9 + n/3. If the maximum value of y; i.e., y mY , is represented 
in the counter by 2 n counts, then 
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e < 


7/9 + a Jl 


,n 


Yrnax 


(2.25) 


Combining Eqs. (2.24) and (2.25) the error bound is 


; 7/ 9 + n/5 


£ (k) ' 2 n W w ( x (k)) %/ x (k) 


- a) 


(2.26) 


This equation has the property that the error bound decrease 
as the number of stages of the BBM is increased. However, it is too 
pessimistic for design purposes. This point will be illustrated 
in the next chapter by applying this formula to a specific countup- 


countdown machine. 
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CHAPTER III 

GENERATION OF FUNCTIONS 


Synthesis (Differential Equation) 

The method of synthesis which will be applied in this section 
is to express the function to be generated as the solution to a 
differential equation* 2t- has been demonstrated with analog tech- 
niques that a wide variety of functions can be generated by 
utilizing only integrating units and adders (e.g. , Ref. 8). For 
example, with a mechanical differential analysis the basic units 
are a ball-disc integrator and a differential. In the synthesis of 
countup-countdown machines, the integrator model of Fig. 2.2 and 
the anti-coincidence circuit which permits the summation of two 
pulse streams will serve as these operational units. It should be 
reemphasized that the principle design elements were recognized as 
entities only for purposes of synthesis, and that the fabrication 
of the actual machines may permit circuit simplification which may 
result in reductions of the hardware requirements. 

The first step in this synthesis is to express the function tc 
be generated as a differential equation such tha ^ the highest order - 
derivative is isolated; i.e.. 



(3.1) 


The independent variable in the above equation is -represented by a 
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clock. The next step in the synthesis is to assume that a circuit 
has been designed to generate the highest order derivative. Inte- 
grators may then be used to successively reduce the order of the 
derivative accoi'ding to the equation 


dk-y 

dx k * a 


<£l d X 


(3.2) 


A circuit for generating the highest order derivative, -whose exis- 
tence had previously been assumed, may now be developed by the con- 
straint defined by the right hand side of Eq. (3.1). The above 
process terminates the design of the basic configuration for gener- 
ating the function. 

The sche ntic diagram of the machine just designed must then 
be scaled in order to (l) exactly match the defining equation and 
(2) accommodate the range of variables in a finite machine. In par- 
ticular, the counters which are used in the machine configuration to 
handle magnitudes having a finite excursion based on the number of 
stages. Therefore, when a bidirectional counter is used then its 
magnitude must be such that 

| counter value j < -2 n - 1 (3.3) 

Since the level setting .of a BRM must be less than one, then when 
a counter is used for this purpose its scale will!'. be ?r educed 
accordingly, i . e . , 

[level setting of BRM] = 2 _n [ counter value] (3.4) 

Finally, the scale of both sides of the defining equation; i.e., 

Eq. (3^1), must be the same. . 
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This procedure may be reduced to a finite number of steps. A 
synopsis of these steps is as follows: 

Step 1 . Isolate the highest order derivative in the differential 
equation as shown in Eq. (3.1). 

Step 2 . Assume the highest order derivative has been generated in 
a counter . 

Step 3 . Generate each successive lower derivative by using an 
integrating unit. 

Step 4 . Constraint the independent variable, the function, and its 
derivatives according to the right hand side of Eq.. (3.1) and con- 
nect its output directly to the counter representing the highest 
order derivative. 

Step 5 . Assign arbitrary constants to the independent variable and 
its highest order derivative. 

Step 6 . Write constraint equations at each counter based on the 
maximum excursion and number of stages. 

Step 7 . Write constraint equation based on the defining equation. 
Step 8 . Calculate scale faevors to satisfy the equations of Steps 6 
and 7. 

We continue with the application of this procedure in the de- 
sign of specific c or uit down - c ountup machines. The first -two exam- 
ples will be a machine for generating the exponential, function, and 

j j 

another machine for generating the sine -cosine functions. These 
two machines will be illustrated in detail. 
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Exponential Function 
The differential equation 

y’ = y, y(o) = l (3.5) 

whose solution is 

y = e x (3.6) 

will be used to design the circuit for generating the exponential 
function. The design solution for this circuit is presented in 
Fig. 3.1. The detail procedure for this synthesis is as follows: 
Step 1 . The defining differential equation given by Eq. (3.5) is in 
the desired form. 

Step 2 . Assume a circuit has been designed to generate the highest 
order derivative.. This is represented by the line labeled y 1 in 
Fig. 3.1a. 

Step 3 . The function y is generated by integrating y ' . 

Step 4 . Since by Eq. (3.5) the assumed highest order" derivative 
y' is equal to y, then y is diiectly connected to- the line y' . . 
This completes the basic circuit shown in Fig. 3 -la.. 

Step 5 . This is the first step in the design of the scaled circuit 
shown in Fig. 3.1b. The arbitrary constants A and B are assigned .as 
scale factors to the - independent variable r and the highest order 
derivative, respectively. The interpretation of A is "A counts 
per unit of x". The interpretation of B is r B counts per unit 
of.- y ,,: Note that the scele factor of the counter - in Fig. 3.. lb.: is 
reduced by 2” n when it is used to set the levels of the BRM. If 
the value of the counter is By’ counts, -then the level setting of . 
the BRM is- 2" n By' . ' 
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Step 6 . Constraint equations are written for the counter; i.e., 

l 2 "“ AB y| aa x S 2“ - 1 (3-7) 

Step 7 . The mechanization of the defining equation is justified by 
the following constraint equation. 

By' = 2" n ABy (3.8) 


Step 8 . Erom Eq. (*-8) it can be calculated that 

A = 2 U counts/unit of x (3.9) 

Bae calculation of B depends on the maximum excursion of the 
variable y according to the equation 

B |ymax| < 2 n - 1 (3.10) 

This essentially completes the schematic design of the circuit for 
generating e x . 

In order to select the number of stages; i.e., n, it would be 
desirable if an equation were available relating n to the accuracy 

•r* 

of the machine. Unfortunately, the equation obtained in CHAPTER II 
is not suitable for this purpose. This may be demonstrated for this 
machine by calculating the bound given by Eq. (2.26) for various 
values of n. Using L = 1 and y^y = e for this process, then : 



e(e - 1) = 1.907 


fl T 5n 
V 2 n 


(3.11) 


Eq. (3.11) may then be used to calculate the bound for this machine. 
These calculated values. are presented in tabular form in Table 3.1 


for various values of n. 



TABLE 3.1 - ERROR BOUND FOR 
EXPONENTIAL MACHINE 


n 

Bound 

6 

1.190 

7 

,892 

8 

.664 

9 

. 49.". 

10 

.362 


This bound does have the property, which was observed earlier, 
of decreasing as n increases. Unfortunately, it does not decrease 
rapidly enough for design purposes. 

Using a value of y^x <3*2, some values of B have been 
calculated and are given in Table 3.2 together with the initial 
value of the counter to match the initial condition y(0) = 1. 

TABLE! 3;2 - SCALE! FACTORS AND INITIAL CONDITIONS 
OF EXPONENTIAL MACHINE 


n 

A 

B 

initial counter 
value 

5 

32 

10 

10 

6 

64 

20 

- 20 

7 

128 

40 

40 

8 

256 

80 

30 


This series of machines have been simulated on a computer and the 
results are presented in Fig. 3.2. When these results are compered 
to the desired output it is immediately observed that these results 
are much better than those predicted by the error hour is given by . 
Eq. (3.11). Moreover, the five and seven stage machines are seen 
to give better results than the six and eight stage machines* 

A more realistic evaluation of the results presented in 
Fig. 3.2 would be to compare them to the difference equation solu- 
tion. The difference equation for the exponential machine may be 




L40 160 


180 


m 


Figure 3.2. - Exponential machine output, 
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obtained directly from Fig. 3.3, which is identical to the circuit 
shown in Fig. 3.1 but labeled according to Euler’s integration. 

The scale factor associated with the counter *3 "B pulses per unit 
of y." When the counter is at iterative step k-1, then its value 
is at y( £_}_)* During this iteration 2 n Ac pulses arrive at the 
BRM counter and the BRM puts out By(it-l) Folses. These output 
pulses are added to the counter to form the counter value for 
iterative step k. Mathematically the value of the counter may be 
expressed by the difference equation 


By (k ) - %( k .i) + By (k . l) Ax (3.12) 

If each clock pulse is taken as an iterative step then Ax = 2 -n 

and Eq. (3.12) may be rewritten as 

y(k> - yoc-u* 1 + z ' n) (3 - 13) 

Solving Eq. (3.13) in terms of the initial conditions of y (i.e.. 


y(0) = l), then 

y (k ) » (1 + 2 _n ) k (3.14) 

The difference between the difference equation solution and 
the differential equation solution ^s the truncation error of the 
process. The difference between the difference equation solution 
and the actual output is error due to round-off and is illustrated 
by the difference in the curves shown in Fig. 3.2. Since the 
round-off error has been shown to be dependent on the starting value 
of the BRM, it can be changed by using a different starting value 
with the objective of obtaining better agreement between the 
solution and the actual output. Fig. 3.4 presents the actual output 
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for the six stage 5BM for a number of starting values together with 
a plot of the desired results. It is noted that by iris simple 
means measurable improvement has been attained in the total error. 

Finally, consider a configuration of this machine built out of 
logic elements cm the Case Logic Breadboard (see Ref. 10) for gener- 
ating the function. This is shown for a four stage system in 
Fig. 3.5. Since the exponential function is a monotonic increasing 
function, the counter show:' in this circuit is a simple forward 
counter . 

The synthesis and analysis of a coun tup-countdown machine for 

the generation of the general exponential function 

y = y 0 < e** (3.15) 

from the differential equation 

y' = a, y(o) = y 0 (3.16) 

follows with only minor modification the design of the machine for 
generating e x . The schematic diagram for this machine is given in 
Fig. 3.6a and the logical design is given in Fig. 3.6b. 

Sine-Cosine Generator 
The differential equation 

y” = -y> y‘(o) = o, y(o) = i (3.17) 

is used to design the sine-cosine generator. The he sic schematic 
diagram and the scaled schematic diagram are shown in Figs. 3. la. 
and 3.7b, respectively, and may be developed systematically as 
follows : 

Step 1 . The defining differential equation given by Eq. (3.17) is 


in the desired form. 










(b) Scaled design. 

Figure 3. 7. - Sin-cos generator. 
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Step 2 . Assume that a circuit has been designed to generate the - 
highest order derivative. This is represented by the line labeled 
y” in Pig. 3.7a. 

Step 5 . The function y' is generated by integrating y", and the 
function y is generated by integrating y’ . 

Step 4 . Since y enters the differential equation negatively, the 
pulses arriving at the y counter have a sign change with the re- 
sult that the output of the counter is -y corresponding to y" 

(see defining equation). 

Step 5 . Arbitrary constants A and B are assigned as scale 
factors to the independent variable and the highest order 
derivative, respectively. The interpretation of A is ”A counts 
per unit of x" (i.e., per radian). The interpretation uf B is 
"B counts per unit of y **. If y’ counter is set to 0 and y 
counter is set to -1 (note that this corresponds to the cosine 
being +l), then y 1 will countup to generate sine and y will 
countdown to generate the cosine. 

Step 6. Constraint equations are written for each counter; i.e., 

AB2-»iy’| mx < 2 n - 1 (3.18) 

A^-^Jyl^ < 2 n - 1 (3.19) 

Step 7 . Constraint equation is written to justify the defining 
equation. 

By” = -A 2 B2~ 2n y (3.20) 


Step 8 . 


Pror Eqs. 


(3.10) to (3.12), 
A = 2 n , Bjyjjjjg^ 


the scale factors can be chosen. 


(3.21) 


< 2 n - 1 
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Sot *» choices for A . and B are given in Tab3,e 3 . 5 cased c . 
j3q. v o.81 j < 

T'lE; r - -.5 - SvA TS S.0TCF6 I1FTIAI. 

ONI’ '’IONS KK sin-cos WKPJS 


n 

a 

r~_ — 
B 

s.n 

counter 

003 j 

counter 

3 

8 

4 

0 

-4 

4 

16 

8 

0 

V. ■ ^ 

5 

32 

16 

0 

-16 

6 

64 

32 

0 

. -32 


This series of machines has been simulated on a computer and the : 
results plotted on Fig. 3.8. ^ , 0 

The truncation error associated with this circuit may be cal- 
culated by solving the difference equations associated with this 
circuit*- Calling the values of the sine counter and cosine counter 
at iterative step k ! ’?£(^ .ana respectively, ' the differ- 

ence equations at these two counters are : 

- 3S (k ) - E3( k-1) +.BC (k .i) Ax ■ / 

-38(^3.) SSik-l).:^ (3.82) 

or mo^e concisely - 

75 )? (? - 23) 

where V^S k is the second backward difference. Therefore, the 
second derivative in this machine is approximated by the second 
backward difference . 

If each cxock pulse is taken as in iterative step, then 
. (3.22) may be written in matrix form as 
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An approximate solution of 3q. (3.24) for lar^e n in terms 
of the initial conditions 



may be written as 

s (k) * (1 + k2' 2n_1 ) sin 2 _n k (3.26) 

C (h) = t 1 + itf' 211 " 1 ) cos 2“ n k (3.27) 

quantitative evaluation of this circuit is complicated by the 
fact that it is used to generate two functions. One method which 
seems especially well suited for testing such a circuit is to plot 
one output function with respect to the other function rather than 
with respect to the independent variable. For the sin-cos generator 
this is called the "circle test" since the resultant figure for a 
perfect sin-cos generator would be a circle. Moreover, it is 
possible to study the errors due to round-off independent of those 
due to truncation by comparing the actuil output to Eq. (3.24) out- 
put. For the sin-cos generator a composite piou of the solution to 
the difference equation may be simply obtained by expressing 
3qs. (3.20) a nd (3.27) in polar coordinates, with the result that 

p = (1 2‘ n-:L e) (3.28) 

The difference between this equation and the circle represents the 
truncation error of the process and is seen to increase as the 
spiral of Archemedes. The results plotted in Fig. 3.8 are compared 



to the plot of Eq. (3.28) in Fig. (3.9) by this method. 

Other Differential Equation Machines 
The sinh and cosh machine is based on the differential equation 

y" = +y (3.29) 

The schematic diagram for this machine is similar to the sine-cosine 
generator except that all outputs from the BRM's are added to the 
counters. The basic circuit and the scaled circuit for this machine 
is shown in Fig. 3.10. It will be noted from this diagram that the 
values of A and B are defined by the equation set (3.30) below. 

By” = A 2 B2" 2n y 
A = 2 n 

^l^maxl — ^ ^ 

B |y-mxl ^ " n - 1 < 3 * 30 > 

The output for a sinh and cosh generator is plotted in Fig. 3.11 
for a five stage system where B is chosen equal to 16. It is 
instructive to display the p-sequenee calculation from which these 
results were obtained. These are shown in Fig. 3.12. 

A series of other useful machines will be illustrated in this 
section. In particular, if two pulse streams du and dv are 
given, then the product, uv, may be generated by using the equation 

duv = u dv + v du (3.31) 

The basic design for this product machine is shown in Fig. 3.13. 

The machine for generating the square of a function is shown 
schematically in Fig. 3.14. This machine will generate the function 
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Figure 3.9. - Sin-cos circle test. 
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Figure 3.12. - P -sequence ciilcu3.ationi fcK sinh-cosh machine. 
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and is based on the equation 

y' = 2x (3.33) 

The square machine is utilized as a subassembly in the machine 
for generating the reciprocal of a function^ i.e., 

y = l/x (3.34) 

The basic design for this reciprocal machine based on the differen-. 
tial equation 

y' - -y 2 (3.35) 

is shown in Fig. 3.15. 

The machine for generating the solution to the second order 
differential 

y" + 2a^,|y’ + a| y = 0 (3.36) 

is shown in Fig. 3.16. 

The tan machine is shown schematically in Fig. 3.17. This 
machine 5s based on the differential equation 

y 1 = 1 + y 2 (3.37) 

A similarity will be noted between this machine and that of the 
reciprocal machine. 

The square root machine is based on the solution of the differ- 
ential equation 

y’ = -l/(2y) (3.38) 

It will be noted by this equation that it will form a subassembly 
of the reciprocal machine. That is, the differential equation 

dz/dy = 2z 2 


(3.39) 
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is used in order to form the function 

z = -l/(2y) (3. 40) 

The basic design of this machine is shown in Fig. 3.18. 

Synthesis (Difference Equation) 

Consider the iterative process of successive substitution in 
the functional equation 

x (k+D = ¥< x (k)> 

where cp(x) is chosen such that the fixed points of cp(x) (i.e , 
the points Xj_ where x i = (pCx^)) are the roots of f(x) = 0. 

One simple form of cp might be x - f(x) which leads to the iter- 
ative process 

x (k+l) = x (k) - f ( x (k)) ( 3 - 42 ) 

A more general form is x - g(x)f(x) which leads to the iterative 

process 

X U+D = *<*) ' g(x (k) )f(x (k) ) (3 - 43) 

A restriction on g(x) in this latter form is that it has no zeros 
that are not zeros of f(x) and that the multiplicity of its poles 
at the zeros of f(x) be less than the multiplicity of the zeros 
of f(x) at these points. With these restrictions, it can be 
readily seen that the iterative Eq. (3.43) has fixed points at the 
zeros of f(x) (i.e., x = x - 0). The function g(x) in Eq. (3.43) 
is chosen so that the process converges. 

The basic equation of a counter immediately suggests a method 
for generating an equation of the form of Eq. (3.43). This method 
of synthesis is simply to generate a pulse stream equal to g(x)f(x) 
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t-ud feed it into a counter. This procedure may be outlined as 
follows : 

Step 1 . Write the defining equation .in the implicit form.g(x)f(x) = 0. 
St ep 2 . Assume a counter with the value of x and generate a 
pulse i_ tream equal to g(x)f(x)p, where g(x)f(x) is the level 
setting of a HRM and p is the input to the BRM counter. 

Step 3 . Feed l-?ck the pulse stream generated in Step 2 into the x 
counter . 

Step 4 . Assign an <- \b? trary constant to each variable represented 
in the machine. 

Step 5 . Write constraint equations and calculate the scale factors 
ouch that these equations are satisfied. 

Divide Algorithm 

The machine for generating x such that 

x - a/b (3.44) 

may be designed as follows: 

Step 1 . ')ne way in which Eq. (3.44) may be rewritten <"ut it 
into the desired form is 

xb - a = 0 ?,3.45) 

Step 2 . Assuming a counter value representing x, a pulse stream 
equal to xb - a Bay be generated. This is shown in Fig. 3.19a. 

Step 3 . The pulse stream generated in Step 2 is fed into the 
counter representing x. 

Step 4 . The schematic of Fig 3.19a is redrawn in Fig. 3.19b, and 
each variable is assigned an arbitrary constant. 


(a) Basic design. 



(b) Scaled design. 

Figure 3. 19. - Divide algorithm. 
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Step b . The constraint equations may be written directly from 
Fig. 3.19b, in which an iteration is taken every clock pulse. 

°l 3 w! s 2 n - 1 

AjSmaxI - 2 n - 1 

®|kffiaxj ^ (3.46) 

Gx (k+1) = Cx U) " BbCx (k )2- 2n + A& 2 " n (3.<t7) 

Suppose for the sake of argument that 

A = B = C = 2 n (3.48) 

Eq. (3.48) implies that 

^ 1 - 2 " n 
a max — ^ ~ 2 n 

b max < 1 - 2" n (3.49) 

and Eq. (3.47) may be rewritten as 

x (k+l) = x (k) - 2-»bx( k ) + a2- n (3.50) 

If Eq. (3.50) converges to a fixed point, x, then 

x = x -2 _n bx + a2~ n (3.51) 

If S( k ) is the iterative truncation error at iterative step k; i.e, 

^(k) ~ ^ ** x (k) (3.52) 

then from Eqs. (3.50) and (3.51) 

£ (k+l) = (1 - t>2~ n )£( k ) (3.53) 

This may be written in terms of &(o) as 

£ (k) = 0- - b2 " n ) k£ (0) (3-54) 

Tliis process will converge if 

lim 6( k ) =0 (3.55) 

k-*oo 

which implies the condition 
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jl - b2~ n | < 1 (3.56) 

for convergence. Therefore, by Eq. (3.56) the process is seen to 
converge. However, if the sign of the outputs of the ERM are re- 
versed such as shown in Fig. 3.20, then an analysis will show that 
the process will not converge. 

Other Difference Equation Machines 
The square root machine; i.e., 

x » Va (3.57) 

may be designed by finding the zeros of the equation 

x 2 - a = 0 (3.58) 

This machine is shown schematically in Fig. 3.21. If the scale 
factor of x and a are both taken as 2 n , then an analysis 
similar to that of the divide algorithm shows the iterative process 

x (k+l)~ °Kk) " 2 ^k) + 2 *** (3.59) 

is generated by the machine. The iteration truncation error for 
this equation may be written as 

e (k+l) = 1 - 2" n (x + x( k )) e( k ) (3.60) 

where x is the solution. A sufficient condition for this process 
to converge is that 

jl - 2“ n (x + x ( k j)| < 1 (3. 61) 

Since 

X jaax <1 (3.62) 

with the scale factor chosen, then the process converges, however. 


had we chosen 



74 



Figure 3. 22. - 
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x (k+l) = x (k) - 2 ~ n ( a “ x (k)) 2 ( 3 - 63 ) 

as the iterative process, then the process would not converge. 

A product machine may also he designed using this method of 
synthesis. In particular, if 

x = ab (3.64) 

then the product may he found by the .iterative process 


x (k+l). ~ x (k) - P< x (k) " ab ^ < 3 - 65 ) 

This machine is shown in Fig. 3.22. If the scale factors for x, 

a, and b are OJ. taken as 2 n , then the machine generates the 
iterative process 


x (k+l) = x (k) • 2 ' n < x (k) - ab > (5-66) 

The iterative truncation error for this machine in terms of 

6(o) is 

«(k) = (1 - 2-»») k S(o) (3.67) 


Synthesis (Regenerative Circuit) 


Consider the schematic diagram shown in Fig. 3.23. The value 
of K is bounded such that 


jKj < 1 - 2~ n 

The output equation for this circuit may be written as 


dz =* K(dx + dz) 

dz K 
dx 1 - K 


(3.68) 


(3.69) 


As K •* 1 in Eq. (3.69) then the ratio dz/dx -* «>. However, 
Eq. (3.68) fixes an upper bound on this ratio such that 


1 - K 


(3.70) 
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Therefore, by using this regenerative circuit, the BRM may act as 
an amplifier. However, if large amplifications are to be considered 
other factors must be taken into account. In particular, suppose a 
ax pulse arrives at the BRM and this generates a dz pulse . The 
dz pulse is delayed by a gated pulse generator and is fed back to 
the BRM. This in turn may generate another dz pulse. This process 
may be continued depending on the value of K. However, each time a 
dz pulse is recirculated the pulse shape deteriorates . For exan pie, 
if leading edge logic is used then the rise time of each pulse will 
be increased until the dz pulse is not sharp enough to be utilized. 
Moreover, enough time must be allowed between the dx pulses to 
permit the maximum number of dz pulses. The maximum value cf K 
usually utilized in these circuits will in general, be less than 
that permitted by Eq. (3.68). 

This point is illustrated by the plot of Eq. (3.69) in 
Fig. 3.24. If -1 < K <; l/2 then at most 1 pulse will be fed- back 
with each input pulse. If l/2 < K < 1 then more than 1 pulse will 
be fed back with each input pulse. Therefore, by fixing the upper 
value of K, the maximum number of feed-back pulses may be re- 
stricted. 

The method of synthesis in this section is similar to that 
used in the synthesis by differential equation. However, in this 
section the differential equation will involve the highest order 
derivative on both sides of the equation; i.e., 



(3.71) 
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In general the equation is written in this form when the highest 
order derivative can not be isolated In particular, if the highest 
order derivative has a non-cons cant coefficient then it may be 
written as Eq. (3.7l) by adding and subtracting a constant from 
this coefficient (see Refs. 6 and 7). The design of the circuit 
based on Eq. (3.7l) implies the use of the regenerative circuit 
since the generation of the highest order derivative involves itself. 

Square Root Machine 

Consider the generation of the square root 

y = -/x (3.72) 

from the equation 

y = 1/2 (3.73) 

The first step of this technique is to write Eq. (5.73) as 

(y - C + e) g = 1/2 (3.74) 

and then isolate the highest order derivative as shown in the 
following equation. 

dy = §1% - (y - c)dy| (3.75) 

The basic and scaled circuit for generating Eq. (3.75) is shown in 
Rigs. 3 . 25a and 3.25b, respectively. Following the same technique 
used to derive a machine based on a differential equation, the con- 
straint equations are written. 

A dy = B dx - A 2 2" n (y - C)dy (3.76) 

is the defining equation for the machine, and 

A|C - < 2“ - 1 (3.77) 


7 $ 




<b) Scaled design. 

Figure 3.25. - Regenerative circuit square root machine. 



figure 3. 26. - Scaled diagram for four-stage square root machine. 
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is the counter limiting equation. From these equations the scale 
factors may be computed by 

A = 2 n /C, B a 2 n " ,;L /C 2 (3.78) 

where C is chosen such that it satisfies the inequality 

1 - 2d < 1 - 2' n (3.79) 

a Jmax 

Based on this design. Fig. 3.26 gives the computed scales and cir- 
cuit for a four stage regenerative circuit square root machine. 

The output of this machine together with that of the desired output 
is shown in Fjg. 3.27. 

Other Regenerative Circuit Machines 
The natural logarithm machine 

y = In x (3. 80) 

may be designed as a regenerative circuit by the equation 

dy = M + (i - dy (3.81) 

The circuit for generating this function is shown in Fig. 3.28. 

The quotient machine 

z = x/y (3.82) 

may be designed as a regenerative circuit by the equation 

“ a O “ a ) dz + z “ dxj (3.83) 

The schematic circuit for this function is shown in Fig. 3.29. 
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Pulses 

Figure 3. 27. - Output of regenerative 
circuit square root machine. 
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CHAPTER IV 

PIECEWISE POLYNOMIAL MACHINES 
Polynomial Machines 

When a function to be generated is available only as empirical 
data (e.g., such as in sampled data systems) then the design must 
be based on generating an approximation function. Various classes 


of approximation functions ^nd techniques for obtaining approx - 
matious have received considerable attention in the literature 
(e.g., see Ref. 14). A particularly convenient form for approxi- 
mating a continuous function is that of a polynomial. The general 
polynomial 


f( x ) = + .Jfor.i xm-1 + 

nx ' ml (m^TJT 


a 2 

. + -X? + a-jx + a 0 (4.1) 


may be generated by the circuit shown in Fig. 4.1. However, it is 
usually the case that all the data is not immediately available for 
generating the function over its entire range, or if it is available 
the polynomial needed to meet the accuracy requirements is of ex- 
cessively high degree. In these applications the requirements of 
the problem may be met by using a series of relatively low degree 
polynomials where each polynomial is used to fit data only in a 
restricted range. Such machines are called piecewise polynomial 


machines . 

Because of their wide spread use in applications (see Refs. 9, 
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11, 15, and 15) and also because they can be realized with relatively 
simple circuits, this chapter will be devoted exclusively to de- 
scribing a series of machines for generating piecewise polynomials 
arising from finite difference techniques. These machines are 
grouped into two broad categories based on applications; i.e., 
interpolation or extrapolation. Each machine of this series will 
generate a low order polynomial fitted to data available at equal 
intervals of the argument. In passing from one segment into the 
next new data is introduced. The form of the data in each case is 
simply generated from the. empirical data. 

In order to facilitate the description of the machines in the 
next two sections, ordinary difference notation will be used. In 
particular, av, is the value of the independent variable where the 
value of the function is obtained; i.e.. 


f(<%.' = f n (4.2) 

The quantity fro represents the spacing of the independent vari- 
able, and are the successive differences which 

may be obtained from lower order differences as follows: 

4:. “ f n+l " f n 


* ‘W. 



A® „ & ■-! . d 1 ' 1 (4.3) 
n n+1 n ' ' 

In formulating the approximation formula which passes through 
the given points, it is converient to display these various differ- 
ences in tabular form as shown in Pig. 4.2. Prom this difference 



Gregory-Newton Backward formula 



Newton-Sterling formula 
Newton -Bessel formula 


Gregory-Newton Forward formula 


(a) Direct path difference formulas. 

A -i A? 2 

Newton -Gauss Backward formula 


Newton-Gauss Forward formula 

A) A?i 

(b) Broken path difference formulas. 



Figure 4. 2. - Difference table and paths of difference formulas. 
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table, alternate forms of the approximation formula may be derived 
dependent on the differences which are utilized; i.e., on the path 
through the difference table. The paths of some of these formulas 
are shown in Fig. 4.2. The form of the formula which will be 
utilized in the following discussion will be such that in each case 
the value of the variable x will vary from 0 to 1 in the interval 
of interest. 

Interpolation 

A piecewise linear interpolator can be obtained by passing a 
linear polynomial through successive pair of points of the function 
to be generated. This scheme is illustrated in Fig. 4.3. A linear 
polynomial is generated which passes through the points Pq and P^. 
At the point P]_, the point P 2 is added to the scheme and a linear 
polynomial is generated which passes through the points P-j_ and P 2 . 
This procedure may be expressed in terms of ordinary differences by 
the Gregory -Newton interpolation formula; i.e. , 

f(<% + x&oo) = f n + £fcX (4 4) 

The first derivative of this formula is 

&aof’(% + x6a>) = V (4.5) 

This interpolation formula may be generated by the linear polynomial 
generator shown in Fig. 4.4a. 

The corresponding equation and its first derivative for the 
next interpolation interval is given by, 

f (°h+l + x6a >) * f n+l + Ai+l x 
Sarf'( %+ i + xSco) = Vi 


(4.6) 

(4.7) 



Figure 4. 3. - Scheme for piecewise 
linear interpolation. 




(b) Function values input data. 

Figure 4. 4. - Machines for piecewise linear interpolation. 





At the end of the first interval of the values of the function 
and its derivative given by equations (4.4) and (4.5) are 

f(<% + Scd) = f n + = f n+l (4.8) 

6<x>f ' (o^ + 5 cd) = (4.9) 

The corresponding values of these two quantities needed at the start 
of the next interval are given by Eas. (4.6) and (4.7) and are: 

f ^ii+l' = f n+l (4.10) 

8oxT (o^+x) = Afc+x (4.11) 

By direct computation it may be verified that in order to proceed 

O 

from one interval to the next, the quantity Z^ (which is Z^ + x - An) 
needs to be added to the setting of the BRM, and the output (i.e., 
the end point of the interpolation interval) need not be modified. 
However, since adders have been excluded as basic design elements, 
the same result may be attained by transferring Z^+x as the set- 

o 

ting of the BRM (since Z^+x = Z^ + z£). Consequently, the circuit 
shown in Pig. 4.4a may be used for piecewise linear interpolation 
of a function by transferring successive first difference as set- 
tings for the BRM in order to proceed from one interval to the next. 

The circuit derived above is well suited for an application in 
which an incremental encoder is used to generate the input data. If 
an absolute encoder is used to generate the primary data, then the 
above circuit may be adapted for this input by using the defining 
equation for first differences; i.e., Eq. (4.3). This circuit is 
shown in Pig. 4.4b. As in the previous case, only one new piece of 
information must be transferred into the circuit in order to pro- 
ceed from one interpolation interval to the next. However, in this 
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case, before the new information, i.e., f n+i >, is transferred as the 
setting of the lower BRM, the value of the lower BRM; i.e., f n+1 , 
must he transferred to the upper one. 

Two piecewise quadratic interpolators will he derived. The 
schem - for the first one which will he called the hack interval 
quadratic interpolator is illustrated in Fig. 4.5. A quadratic 
polynomial is generated which passes through points Pq, P^, and 
Pg. This polynomial is used to generate the curve between points 
Pq and P^* The point P 3 is then added to the scheme and a 
quadratic polynomial is derived which passes through the points 
Pj_, P2, and P3 This polynomial is then used to generate the curve 
between P^ and P2 . 

This procedure may also he conveniently expressed by the 
Gregory-Newton quadratic interpolation formula ; i.e., 

f(<% + x 8 o>) = f n + ^qX + — X - ~ ■ 4 (4.12) 

The successive derivatives for this formula are 

Sojf ' ((% + x&d) = ~ 4) + 4 X (4.13) 

(8a>) 2 f"(a^ 1 + x8cd) = 4 (4.14) 

The corresponding equation and its derivatives for the next 
interpolation interval are: 

f(oh+i + x&o) « f n+1 + - | 4+1 ) x + x2 U-15) 

8 rf'((%+i + x5o>) = ^ +1 - | 4+1 + 4+l x (4.16) 

(Sa)) 2 f"(a^+x + x8a>) = 4+1 


(4.17) 




Figure 4. 5. - Scheme for piecewise qua- 
dratic Interpolation (back interval). 
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Consequently, the correction to be added to the second deriva- 
tive, first derivative, and function at the end of i;he interval in 
order to proceed to the next interval are -i A^ + ^, anc! 0, 

respectively. From this formulation, however, addition is required 
in order to proceed from one interval into the next. A formulation 
of this process which leads to the elimination of the explicit 
adder is to splinter the polynomial given by Eq. (4.12) into the 


following two polynomials. 

/ V* 

f a (% t x8cd) = f n .+ Ai x + ~2" xZ (4.18) 

/ Vjj 

f b (o^ + x5co) = X (4.19) 

where 

f(o^ + x8o>) = f a (d^ + x8cd) + f b (d' n + x8cd) (4.20) 

The first and second derivatives of Eq. (4.18) are 

(boiJfa'Ca^ + x&cd) = + A^x (4.2l) 

(&co) 2 f_(<D +:xScd).= A? (4.22) 

» n B 

The corresponding splintering of Eq. (4.15) yields 

aLi 

f aK+l + x8c °) " f n+l + Vl x+ ' - T“ x2 < 4 * 23 ) 

(5cp)fa(%+i + x8c °) = Ai+1 + 4+l x (4.24) 

(So>) 2 f a (%+l + x8co) = Ag +1 (4.25) 


Consequently, no correction need be adder to the first deriva- 
tive in generating the function f a . Eqs. (4.18) and (4.19) may 
then be used to design the circuit shown in Fig. 4.6a. It will be 
noted that in order to proceed from one interval to the next, only 
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new second difference date, need be transferred to set the levels of 
the leftmost BRM. . 

Based on the defining equation for second differences; i.e., 

Eq. (4.3), the machines shown in Figs. 4.6b and 4.6c may be obtained 
directly from the machine shown in Fig. 4.6a. In these machines, 
previous first differences and values of the function are trans- 
ferred directly from the lower BRM’s to the upper ones before a 
new first difference and a value of the function, respectively, is 
transferred into the -lower- one in order to proceed from one inter- 
val into the next. 

The scheme for piecewise front interval quadratic interpolation 
illustrated in Fig. 4.7 may be derived by use of the Newton-Gauss 
interpolation formula given in the following equation. 

f (<% + x&d) = f n + x ^ 4-1 U-26) 

If Eq. (4.26) is implemented directly, then the first deriva- 
tive and second derivative must be corrected by adding — A^ 
and A^, respectively, to these quantities in order to proceed 
from one interval to the next. The explicit need for an adder may 
be avoided in a manner similar to that used in the previous discus- 
sion by splintering Eq. (4.26) into the following pair of equations. 

f a(% + = f n + x + | x2 4-1 (4.27) 

f b K + xScd) = | 4-1 x (4.28) 

Based on this pair of equations, the circuit shown in Fig. 4.8a 
may be derived directly. The machines shown in Figs. 4.8b and 4.8c 
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X 

Figure 4. 7. - Scheme for piecewise qua- 
dratic interpolation (front interval). 



















are adapter from the circuit shown in Pig. 4.8a by using the defini- 
tion of the second difference given in Eq. (4.3). 

A cubic interpolator may be obtained from the Newton-Gauss 
interpolation; i.e. , 

f(c^ + x6cd) = f n + x i x(x- l) 4 - 1 + x(x - l)(x +l) ^n_i (4.29) 

However, in this case the discussion will be limited to using this 
formula for central interval interpolation only. This scheme is 
illustrated in Fig. 4.9. The points p 0 > P l> P 2 -» and P 3 are used 
to generate an interpolation formula for interpolating between P^ 
and Pg • The point P ^ is then added to the scheme and the points 
? 1 , P 2 , P 3 , and P_j. are used to interpolate between points Pg 
and P 3 . 

Eq. (4.29) may be applied directly to yield a central interval 
cubic interpolator. However, in this case the third, second, and 
first derivatives must be corrected by adding A 4 ^, 0 , and i A 4 ^ 
to these quantities, respectively, in order to proceed from one 
interval to she next. 

A configuration may be obtained which conforms with the design 
practice of not using an adder by splintering Eq. (4.29) into the 
following pair of equations. 

p a(°h + X ^ CD ) = f n + x |Ai " 2 ^h- 1 ^ + If* ^-1 + ~ ^-1 ( 4 - 3 °) 

f b (co n + x 8 o>) = -| A ^-_ 1 (4.31) 

Eased on this pair of equations, a circuit may be obtained such 
that the function, its first derivative, and its second derivative 
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Figure 4.9. - Scheme for piecewise 
cubic interpolation (central in- 
terval). 
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need not be changed in order to proceed from one interval to the 
next This circuit is shown in Fig. 4.10a. The circuits presented 
in Figs. 4.10b, 4.10c, and i.lOd are modifications of this circuit 
based on the definition of the third difference. 

Extrapolation 

Extrapolation presents an added problem in that the output . f 
the machine must also be corrected in order to proceed from one 
interval to the next. This is illustrated in Fig. 4.11 for linear 
extrapolation. A linear polynomial through Pq and Pp is used 
to extrapolated the values from Fp to P 2 • The point P 2 is 
then added to the scheme, and a linear polynomial through Pp and 
?2 is used to extrapolate che next interval. The predicted value 
P 2 said the new value Fg can be expected to be different. Con- 
sequently, the output must be corrected for this new value P 2 . la 
order to avoid putting a jump in the output function at this point, 
the scheme which will be employed is to put the correction in 
linearly over the entire next interval. This scheme (as well, as 
that of quadratic extrapolation which will be described next) is 
closely related to the Porter -Stoneman digital filters (see 
Ref. 13) and may be extended accordingly. 

The Gregory-Newton backward finite difference formula may be 
used to design the linear extrapolation machine; i.e., 

f(<% + x5cd) = f n + ^_pX (4.32) 

The corresponding formula for extrapolating the next interval is 

f(%+l + x6o>) = f n+1 + Z^x (4.33) 















100 



(d) Function values Input data. 

Figure 4. 10. - Continued. Machines tor plecawlse cubic Interpolation (central interval). 
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X 


Figure 4. 11. - Scheme for piece- 
wise linear extrapolation. 
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Zi these formulas are applied directly , then the circuit must 

O 

be corrected at the end of the interval by adding to both the 

function and its first derivative in order to proceed into the next 
interval. This, however, would cause a jump in the output function. 
This jump may be avoided by putting the correction term in the out- 
put in a linear manner over the entire next interval. The result- 
ing polynomial has the property that its initial value corresponds 
to the end point of Eq. (4.32) and its final value corresponds to 
the end point of Eq. (4.33). A polynomial which satisfied these 
constraints may be written as follows: 

f ('%+l + x6c °) = ( f n + ^n-l) J \i + 4-l) x (4.34) 

The second difference in Eq. (4.34) may be eljsn ted by 
using the defining equation given by Eq . (4.3). This substitution 
yields the following equivalent equation. 

f(%+l + x&o) = (f n + Vl) + ( 2 Ai ' Vi) x U.35) 

Eq. (4.35) may be implemented to yield the linear extrapolator 
shown in Fig . 4.12. 

The scheme for quadratic extrapolation is shown in Fig. 4.i3. 
The quadratic equation through points Pq, P^, and Pg is used to 
extrapolate the data to the point Pj. The point P 3 is then 
added to the scheme and can in general be expected to be different 
from P* . The quadratic equation through the points P-^, Pg, and 
P^ is then used to extrapolate to the point P 4 . In order to 
avoid putting a jump in the output when new information is added 
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Figure 4. 13. - Scheme for piecewise qua- 
dratic extrapolation. 
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to the scheme, the correction may be put in the output a 
linear manner over the entire next interval in the same manner as 
that employed for linear extrapolation. 

The Gregory-Newton backward difference formula forms the 
basis of the quadratic extrapolation. This formula may be written 
as follows: 


f(c% + x&o) = f n + x V-! + ^ 4-2 


(4.36) 


The corresponding formula for the next interval is: 

ftco^ + x5o>) = f n+1 + x ^ + ZULlJl (4.37) 


If Eqs. (4.36) and (4.37) are implemented directly then the 
quantities 4-2 > 34 _g/2, and 4-2 ^st be added to the output 
function, its derivative, and its second derivative in order to 
proceed from one interval to the other. As was indicated earlier, 
the jump in the output function can be avoided by putting in the 
correction over the entire next interval. A polynomial which 
satisfies these constraints (i.e., has the end of Eq. (4.36) as 
its initial point and the end of Eq. (4.37) as its final point) 
may be written as follows: 


f(%+l + xOco) = f n + -Vl + 4-2 + 





x 


+ 



X 


(4.3b) 

Substituting the difference relationship given by Eq. (4.3) 


into Eq. (4.38) yiel 
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+ x&a) ) = f n + An-1 ^ 4-2 

- (vi - 1 4-i - 4U)* + % x 2 

Eq. (4.39) may be splintered into the two equation 

4 1 

^a^^+l + = + Aa-i + Aa-2 + Ai-l x + 2 ~~ x 

f bK+l + x8a >> “ (| 4-1 - 4-2) x 

to yield the circuit shown in Fig. 4.14. 

The circuits shown in Figs. 4.12 and 4.14 may be readily 


(4.39) 

(4.40) 

(4.41) 

ex- 


panded by use of finite difference relations (as was done for 
interpolators) to yield circuits which accept . *»txonal values 
and first differences as the primary source of data. 
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CHAPTER V 


CONCUJSION 

Summary 

A handful of circuits have been reported in the literature 
which have been designed to meet the needs of special purpose 
digital computer problems arising from real time applications. 

The organization df these circuits is to utilize simple counting 
techniques as the basis for computing and result in a simplicity 
of hardware which make them attractive for such special purpose 
applications. The design of these circuits have been examined in 
this study with the objective of (l) "explaining" the circuits and 
(2) generalizing the design philosophy such that new circuits may be 
admitted with the same organization. In order to be specific 
we limited the principle design elements to three fundamental 
units. The elements are (l) the binary rate multiplier which is a 
means of s;' ling down a pulse stream to some specified fraction, 
x ) oho counter, and (3) the anti-coincidence circuit which is a 
means of separating pulses arriving at a counter simultaneously . 

These design elements are represented as operational units 
which may be used to describe the machines. Operational techniques 
are then used as the method of synthesis. In particular, a counter 
is utilized to represent a first order difference equation and a 
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counter in cascade with a BRM is utilized to represent approximate 
integration. The computational errors; i.e., rounding-off error 
and truncation error, introduced into the machines as a result of 
treating the principle design elements as operational units are 
identified and studied in detail. The rule cf round-off, which is 
simply stated in conventional computers, is not as easily formu- 
lated in these machines. Definite results, however, were obtained 
and the rounding-off error was shown to be dependent on the start- 
ing value of the BRM counter as well as n, the number of stages. 
The approximate error bound of 7/9 + n/3 for the generated 
round-off error proved to be disappointingly pessimistic for pre- 
dicting the propagated error for design purposes. Nevertheless, 
having identified these two sources are errors permitted us to 
obtain better results experimentally by two methods; (l) in- 
creasing the number of stages and (2) changing the round-off error 
by changing the starting value of the BRM counter. 

The method of synthesis is presented in three parts; (l) ex- 
pressing the function to be generated as a differential equation, 
(2) expressing it as the fixed point of an iterative process, 
and (3) expressing it in terms of a regenerative circuit which 
is presented. The method of synthesis is explicitly stated and 
is satisfactory in that all known cire .its may be directly obtained 
from it. A wide variety of other iVr.ctj.ons are also obtained using 
these synthesis techniques. Many of these examples are illustrated 
and in some cases actual experimental results were obtained and 
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discussed with the machine. 

A series of machines is presented for interpolation and ex- 
trapolation of a function which is available only as empirical 
data. In particular, the function is generated over its entire 
range by a sequence of low order polynomials. Finite difference 
techniquss are used to describe the polynomials. The order of 
the polynomial is limited to a cubic for interpolation, and a 
quadratic for extrapolation since these seem to be the important 
cases in practice. Nevertheless, these techniques can be easily 
extended to Include higher order polynomials. 

Recommendations for Further Investigations 

We feed that the choice of prlncip e design elements has 
been correctly limited to units that operate as incremental devices . 
It would be interesting to investigate other components in this 
framework. In selecting the new components two approaches appear 
apparent. First, the components used in this study may be sub- 
divided into smaller functional units with a view of studying sim- 
plification methods of the final design. Secondly, new functional 
units may be introduced with a view of admitting new machines. 
However, if components which operate on the whole word are included 
they should be simple decision type circuits (e.g., sign and mag- 
nitude comparators ) and not new fundamental units like an adder 
which would dominate or supplant the other components. 

It is expected that using worst case conditions for obtaining 
error bounds would not produce satisfactory design results which 
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may be used in the -whole spectrum of problems. The study which 
should produce good results would be to consider each circuit 
individually ana obtain deterministic design results which can be 
applied to that circuit. In particular, an algebraic approach as 
was used in Appendices A and B might yield satisfactory results 
for processes which involve only addition and multiplication such 
as the generation of polynomials. For transcendental functions 
the bounds may be obtained by experimental techniques or perhaps 
by comparing the desired function to one which is attainable by 
algebraic means. 

New methods of synthesis should also be sought either to 
include the pathological cases discussed earlier in the report 
or to exclude them as possible machines. 

Other piecewise curve fitting machines should also be stuided, 
especially those in vhich functions other than low order polyno- 
mials are used and those in which the first and higher order 
derivatives are kept continuous . 

The investigation of these machines would be facilitated if 
hardware and good display facilities were available which would 
permit the circuit to be easily fabricated and studied. We do 
not have in mind the design of still another general purpose com- 
puter since it is felt that these circuits best serve the needs 
of special purpose applications. 
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APPENDIX A 

MULTIPLICATION ERROR BOUNDS (ZERO STARTING) 

In this section the error equation of an n-stage BRM whose 
counter starts with zero will he analyzed with the objective of ob- 
taining tight error bounds. Nevertheless, some of the intermediate 
results which will be obtained in this section are interesting in 
their own right. Because this analysis is involved, we will proceed 
formally. The basic outline is to use Eq. (2.8) to find the points 
where the ir ximum positive value is attained and then evaluate the 
equation at these points. In a similar nr 'ner, Eq. (2,10) will be 
used to find the minimum negative value. 


We begin by stating and proving Lemma A.l. 

Lemma A.l A sufficient condition for E given by Eq. (2.8) to 
attain its maximum, value is that x^ = y^. 

Eq. (2.8) may be rewritten as a bilinear expression such that 
the terms which are dependent cn either x^ or y_^ are grouped 
together. The quanity A in the resultant expression is independ 


ent of either x^ or y_^. 


E(Xi,y_i) - A + 1 Xi y_ ± - | xi^| y.i»! + ~ y„i»2 + • * • + y- n ) 

1 (l ^ 1 . . 1 \ /. ..v 

' 2 y -i(2 x i-l + 4 x i-2 + ' ' • + ilJ X lJ <A ’ 1) 


(A.l) 


By direct evaluation the value of this expression is as follows: 
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E(0,0) = A 

E(0,1) = A - I (I y-i -1 + • • • + y-n) 

E(1,0) =A-|^x,. 1+ • ■ • + xjj 


Moreover, for the specific case when i is n then the value of 
Eq. (A.l) is 


E(0,0) = A 
E(0,l) = A 

E(1,0) = A “ | {j x n -l + * 
E(l,i)=A+ + 



This lemma is proved by observing that the value of E when 
x^ = y_^ = 0 is always greater than or equal to the value of E 
in both cases when ^ y_^ and that for the nth component the 
value of E is always greater when x n = y_ n = 1. 

Based on this lemma, the maximum value of Eq. (2.8) will be 
found by finding the maximum of the quadratic form expression 


QCx^Xg, . . .,:x^) = (x^Xg, 



Theorem A.l For all values of the components x^. 


(A.2) 


Q(1,x 2 ,x 3 , . . .,3%) > Q(0,x 2 ,x 3 , . . ^Xq) 


Eq. (A. 2) may be rewritten such that A is a quadratic ex—.. 


press ion independent of xj_. 



lie 


^ X 1 X 

QCx^Xg, . . .,x n ) * A + x^g- " 4 x 2 " 8 X 3 * * * “ 


2 n 


*n 


(A. 3) 


By direction evaluation of Eq. (A. 3) 

Q(0,x 2 , . . .jXn) = A 

Q(l> x 2> • • *> x n' = A + ^2* ~ 


1 

4 x 2 




Since 


11 1 
2 “ 4 x 2 “ q x 3 


an 


x n > 0 


then 


Q(l, Xl ,x 2 . 


^Xjj) > Q( 0, Xl ,x 2 , . . 


theorem A, 2 For all values of the components 

Q(x^>x 2 , . . ♦ ' > Q( x i* x 2 > * * *> x n-l , °^ 

This proof proceeds similar to Theorem A.l. Q may be re- 
written as 

Q(x 2 >x 2 , . . . >Xn )=A .+. - 4 x n-l “ * • • “ j^n x ^ (A. 4) 

where A is independent of x n . 

By direct evaluation 

Qf x l' x 2 > . . .,x n _ 1 ,0) = A 

Q( x 1, x 2? . . * ,x n- 1 , l) = A + 2 “ 4 *n-l “ * * * “ ^n X 1 
Therefore 

* * * ,X n-l ,:L ^ > * * * ,X n-l ,C) ^ 

Theorem A. 3 For all values of the components x i# 
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Expanding Eq. (A. 5) and. using Eq. (A. 7) to simplify the cross 

product terms, then a typical term x^ may be written as 

( 1 1 1 1 1 . 1 \ 

x i\2 ” 2 i " 2 n-i+l " X m 8 " * * * g?.- i J . 


- x i 


= Xj 




1 _ _ 1 . 1.1 

gi 2 n-i+l 4 8 


2 1 


*) 


f + x i ) |i54fl '|rj 


(A.8) 


Since (xj^ + x^) = 1, then Eq. {A. 6) is independent of the x ± 
variable and the contribution from this term is 


(A.a) 


2 n-i+l 2 i 

Similarly, the contribution to the difference expressed- in 0 

, / 

Eq. J[A.S) from the term invol' '‘r<g x n _^ + ^ is - 

- ' \,“(A.10) 


1 


u 


k . ... i - > 

,1 gn-i+l 


j 


Consequently, the contribution to the difference expressed by 
Eq. (A. 6) by each element may be 'paired by the contribution from 
another element such as to cancel each other out of the expression. 
If ;n . is. odd, then +he middle term .cannot be paired.. But! since- 
this is the (n + l)/2 term, then by Eq. (A. 9) its contribution is . 


1 , 


1 


2 (»+D/2 ~ 2 (n+l)/? f * S ° . ' 

- .j ^ / • 

Therefore, the value of the difference shown by'Ea. (A.S). is ... 
equal to. zero. This implies that - 

' > - - U 'V i‘, ' ' 

' ■ , - ' il J _ ■ 

Q(i,Xg, • •- 


where v and X are the cocrponem by component complement of v 
end a, respectively. 

By Theorer. A. 5 it is noted that 

(Therefore, the difference between the tvo quadratic forms of 
the lemma can be -^pressed as 

6 = Q(l,v,0,a,l) - Q(i,v,0,a,l) (A. 13.) 

Partitioning the M matrix of Eq. (A.U) such that and 

Mg are compatible with the vectors, then 5 may be written as 


5 = (l,v,0, a> l)yYvV (l,v,0,a,l)^^) 


~ (I^O^iJmJV ) - (l,v,0,a,l)%(^ ) 

W U 


(A. 12) 


But it will be noted that 


(l,v,0,a,l)M 1 [v r ]= (l,v,0,a,l)M-j( V 

\0j 


Therefore 


» (l,v,0,a,l)l%^ - (l,v,0,a,l)%^ (A.13) 


It will now be proved by induction on the length of a that 
6 > 0. It may be immediately verified that 8=0 for a of 
length zero. Assume that 6^ > 0 where 8^ is the value of 8 
when a is of length 2k. It will now be verified that 8^^ <1 0. 
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/ x 5 

• I 


f . 1 

• 1 

• 1 

xi 1 

Oi+l 1 

• 1 


Oi+2 1 
li+3 1 

• 1 

* 

f 

• f 

♦ I 

1 i+2k-l 



Oi+2k 

* 


li+2k+lj 

• 

1 

l 

Oi+2k+2l 



1 i+2k+3/ 

U/4 / 


1 i+2k+3/ 

U/2 / 

\ 


1 1 /_i/ 2 l+2k+2| 

x 2 | / -3/2 1+2K ' i ‘ 1 


xi ! : -l/2 2k+3 
Ol+l i : -l/2 2k+2 

o i+2 ; -i/2 2k+1 

li+3 ! -l/ 2 ^ k 


Oi+2k 

1 i+2k+l 

li+2k+2 



-l/8 
-1/4 
| 1/2 


* 3 I 

1 

1 


/ 


(A. 14) 
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The O's and I's in Eq. (A. 14) are subscripted to show their 
position in the vector, and the vector itself is displayed as a 
column rather than a row in order to chow the correspondence be- 
tween the terms which must be multiplied to form the value. 

Multiplying out the terms of Eq. (A. 14), then 

1 


"k+i 


1 1 L 1 1 \ 

Tc 2 ~ 8 ^ 4 * * * ' 2 2k J " 


1 (r 


?2k+5 


2 i+2k+4 

(** + | t-l + - • - + 5b) 


■ 1 
2 2kH 


3 + -L 

4 16 


f + - • • 

+ i 

2 2k-2 

1 tr 

• 2 x i-l ■* 

• • + 


? i+2k+2 

a 

_\ 


(l + J- + . . .+ y ■ ) - — - 

\ 4 J 


i+2k+4 


2 2k+S 

Using the relations 


1 + ~ + 
4 


f/ X jl \ 

^ X i + 2 X i”l + • * * + ~TZ2 X2 J 


• JL - if’ 

• * • ^ 2 2k - 3 [ x - 




and 


1 + 4 * * ‘ * + 


L_ . i | x . JL\ 

l- 2 3 ^ 2 2kJ 


then Eq. (A. 15) becomes 


3k+l = 5 k + 


i-2 

1 .V X 


2 2k+3 2 i+2k+3 2 2k-r4 / , 2 j X i-J 

j=0 


(A. 15) 
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i?>t by direct evaluation 



, _1 1^ 2 1 

2 2k+4 2 i-l 2 2k+4 ~ 2 2k+3 


Therefore 8. + , > 0. 


Lemma A.3. For v = xg,X 3 , . . .,xj- and a = 1, 0,1,0, . . .,1,0 

Q(l,v,l,0,a,l) > Q(l,v,0,0,a,l) 
where v and a are defined as in Lemma A. 2. 

Proceeding in a manner similar to Lemma A.2, it is first noted 
by Theorem A. 3 that 

Q(l,v,0,0,a,l) = Q(l,v,l>l,a,l) 


Therefore, the difference between the two quadracic forms of 
the lemma can be expressed as 

8 = Q(l,v,l,0,a,l) - Q(l,v,l,l,a,l) (A. 16) 

Partitioning the M matrix of Eq. (A. 16) such that and 

M 2 are compatible with the vectors, then 8 may be written as 

a 


A\ , o\ 

= (l,v,l,0,a y l)lw v. J + (l,Y,l,0,a>l)}y a! 

W W 


- (l,v,l,l,a,l)Mjl?J - (l,*v,l,l,a,l)M 2 f a 


(A. 17) 


But it will be noted that 


Therefore 


(l,v,l,0,^,l)^^(l,v,l,l,a,l)M^v 


8 » (l,v,l,0,a,l)iy aj - (l,v,l,l,a,l)}Ja 


(A. 18) 
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Using 5^ to denote the value of 6 when a is of length 
2k, ic will now he proved by induction on the length of a that 

^ x i + g2k+4 x i-l + ♦ ♦ * + 7i+2k+l x 2 * 2 i+2k+2 > 0 


First we note that for a of length zero 


(A. 19) 


8 0 ” 


/l \ 

!-l/2 1+3 \ 


I 1 \ 


/- l/2 i+2, 


1 1 

-l/2 i+3 

! *2 j 


1 


X 2 


[• 


- z 

~l/2 i+2 

5c 3 1 

• 


• 

. 


x 3 

• 




x 3 

♦ 

• 

• 

• 

• 

Xi 


♦ 

-1/16 


• 

% i 


-I/ 8 i 


* 

^i 

• 

-1/16 

lx+1 ] 


-1/8 


i 1 i+l j 


-1 A 


. li+1 

-1/8 

0i+2 


-l/ 4 1 


li+2 


1/2 


li+2 j 

-l/4 

li+3 j 


l/2 1 


\li+3/ 


\0 / 

\li+3/ 

\ 1/2 


p- xi + P X i-1 + 


+ X — XO + , r-T-r— 

2 i+ l 2 1+2 


(A. 20) 

The notation used in Eq. (A. 20) is similar to that used in 


proving Lemma A. 2. 

Assume that 6^ > 0, it will now be shown that 5^+1 > 0 
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8 k+l = 8 k+ 



1 i+2k+l 
0i+2k+2 
li+2k+3; 
\ °i+2k+4| 
\li+2k+g 


_l/2.i+2k+5« 

-i/2 i+2k+4 i 


-l/2 2k+6 


' _i/gi+2k+4> 


-l/8 
. -1/4 
\ 1/2 


* 


;l ■ 

> 1 ' 

L * . 

) • -1/4 

\ J- i 1/2 

\ 1 i+2k+5 \ 0 


Jlj -l/2 i+2k+5 J 

I ~ 2 i I 

x 3i ; 



j 1 /2±+2k+% 

!' 



+ 


Xi 

li+1 

li+2 

Oi+3 

ii+4 


* J i * 

[ Oif2k+l,' 1 “• 
\li+2k+2j \-l/4 
\H+2k+^' \ l/2 


(A. 21) 


Multiplying out Eq. (A.2l), then 0 ^+^ becomes 
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S k+1 = 8 k + 


e 1 + ^ + 


3. 


1 - 


^ “ o?V+R x i 


2 2k+2 i 2 2k+6 


2 i+2k+4 


x 2 - 


1 3 3 

2 i+2k+5 ”4 16 


i + i+ +JL 

i ' 22k 


)- 


2 2k+e 


* plf6 ^ + 


3 - . 3 

x 2 + 


2 i+2k+4 


*1 


2 i+2k+5 2 


| (l + 4 + • * • + 

4 \ 2 “ 2 


! \ 1 4 _ 

2k J 2 2k+3 “ 2 2k+6' ** " 


x 2 - 


2 i+2k+4 “■* 2 H 2k+5 

Equation (A. 22) may be reduced to 

Sk+l = ^ “ ,2k+5 x i ” ’ * 


(A.22) 


1 - 1 

>i+2k+3 X 2 r 


2 i+2k+4 


(A. 23) 


Using Eq. (A. 19), then the right hand side of Eq. (A. 23) 
becomes 

1 — 1—1 1 — 

> 2 2k+3 x i + • • • + gi+gjc+x x 2 + 2 i+2k+2 ~ pk+5 x i “ * • • 


1 - 1^1-, 

x 2 “ „^4-9V4 ^i x i + • • • 


1 *2 + . 1 , ., > 0 


2 i+2k+3 * 2 i+2k+4 2 2k+5 ’ * * ’ 2 i+2k+3 ' 2 i+2k+4 

Theorem A. 4. There exists a v* such that for all v 

Q(l,v*,0,a,l) > Q(l,v,x 1+1 ,a,l) 

where v and a are defined as before. 

First, we note that Q(l,v**,0,a,l) > Q(l,v,l,a,l) because sup- 

* 

pose it were false then there would exist a v** such that for all 


Q(l,v**,l,a,l) > Q(l,v,0,a,l) 


of v 
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But from Lemma A. 2 

* * 

Q(l,v**,0,a,l) > Q(l,v**,l,a,l) 

Therefore, a contradiction exists. 

Moreover, we note that 

'ft# 

Q(l,v**,0,a,l) > Q(l,v,0,a,l) 
that is, there is a largest. 

Therefore, Theorem A-,4 is proved by choosing either v** of 

** 

v** for V*; that is, whichever make Q(l,v*,0,a,l) the largest.. 
Theorem A. 5. There exists a V* such that for all v 

Q(l,v*,l,0,a,l) > Q(l,v,xj +1 ,0,a,l) 

First, we note that Q(l,v**,l,0,a^l) > ‘Q' v l,v,0,0, a,l) because 
suppose it were false then there would exist a v$* such that for 
all v 

Q(l,v**,0,0,a,l) > Q(l,v,l,0,a,l) 

But from Lemma A. 3 

Q(l,?**,l,0,a,l) > Q(l,v**,0, 0,a,l) 

Therefore, a contradiction exist. 

** 

Moreover, a v** can be ehooscn such that 

"ft# * 

Q(l,v**,0,C,a,l) > Q(l,v,0, 0, a,l) 

Therefore 

Q(l,v*,l,0,a,l) ^ Q(l,v,x^ + -^,0,a,l) 

Theorem A. 6. There exists a v* such that for all of - v 

Q(l,v*,0,l) > Q(l,v, 

x i+l>l) 

first., we note that there exist a v** such that 
Q(l.v**,0,l) > Q(l,v,l,l), because suppose it were false then there 
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would exist a v^* such that for all v 

Q(I,t**, 1,1) > Q(l,v,0,l) 

But by Theorem A. 3 

Q(l,V**,0,l) « Qj[l,V**,l,l) 

Therefore, a contradiction exists. 

Moreover, a v^* can be chosen such that 
Q(l,V**,0,l) > Q(l,v,0,l) 

Therefore 

Q(l,v*,0,l) > Q(l,v,x i+1 ,l) 

r w > 

It will now be demonstrated by an example how these theorems , t 
can be used to obtain the Value of x such tha-p the error is the 
maximum positive value. Suppose we consider a 7 -stage EPM. By 
Theorem A.l and Theorem A. 2 

/ 

Q(l,x 2 ,X3,X4,x 5 ,X6,l) > Q(x 1 ,x 2 ,X3',x 4 ,X 5 ,x 6 ,X7) - 

by Theorem A.6 - . - . • 

Q(l,3C 2 ,.3^,X*,Xg, 0,l) ^ Q(l,X 2 ,X2,X 4 ,Xg,X0,l) ' 

by Theorem A.5 

by Theorem A.4 

Q(l,x 2 ,X3 ,0,1,0,1) j ' <qi(l,-x 2 *,X3 ,x^* ,1,0,1) 
by Theorem A.5 - 

Q(l,x§*,l,p,l,0,i)-> &(l,xf*,x§*,0, 1,0,1) 

Theorem A.6 

** 

Q(l, 0,1, 0,1, 0,1) > Q(1 ,x 2 ,1,0,1, 0,1) 
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Putting these inequalities' together 

Q(l, 0,1, 0,1, 0,1) > Q(x 1 ,x 2 ,x 3 ,x 4 ,x 5 ,x 6 ,x 7 ) 

Moreover, "by Theorem A. 3 

Q(l,l,0,l,0,l,l) > Qj[x 1 ,x 2 ,x 3 ,x 4 ,x s ,y 3 ,x 7 ) 

Using the theorems in the pattern illustr- ed by uhe example, 
it is easy to verify that the maximum positive value will occur at 
the points shown in Table 2.4. 

The positive value of .the error' may be expressed con- 

cisely as follows: 

Let denote this -.value for an k. stage EBM* If. } 


*5 . 


• f. 




is odd, then 




f-1 1 

-i/2 k+ 2] 

02 ? 

. 1 

13 ' 

• • 

; * ' . 

lk 

-1/8 : 

Qfc+1 

- 1 / 4 ' 

lk+2 

l/2 


. (A.24) 


Evaluating Eq.. (A.25) J , this yields the difference equation 

*«<**?) •- ye* 

. . \ Ci • _ • - v , 

Solving Eq.. (A.30) for a n stage ERM in -terms of 
yields: 


' / 
/ 


' £ * 6 ” Q.pn 


a* W l) = 1/3 
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Therefore, the maximum positive errer of an n-stage BRM 
where n is odd is: 


E 

max 


(a) 



\A.25) 


If k is even, then 




i Wx( k + 2 ) = + 


03 

14 


°k-l 

Ik 
• °k+l 

\ lk+2 i 


- 1/8 

-1/4 

1/2 




1 1 
3 gk+2 


(A.27) 


Solving this difference equation for an n-stage HIM in terms 
of ^, wy (2) and then evaliiating the r,sultant expression for 
3««(2) - 3/4, yields 


- nr + 1 + ~n (A - i:8) 

Combining equations (A.26) and (A. 28) gives a closed form 
equation for the maximum positive error of a n-stage BRM. 

*kax( n ) = 18 + 6 + gV 2 n (A. 23) 

By applying Eq. (2.10) the minimum negative value for a n 


stage BRM can be obtained. Coinparing the form of Eq. (2.8) to 
Eq. (2.10), it is seen that our previous results can be utilized 
with a slight modification. In particular, the value of the minimum 
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is equal to the negative of the maximum and occur at points which 
are the 2 *s complement of the maximum value. Consequently, the 
minimum negative values will occur at the points shown in Table 2.5. 
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APPENDIX B 

MULTIPLICATION ERROR BOUNDS (ARBITRARY STARTING) 

The error formulas given by Eqs.(2.12) and (2.13) expresses 
the multiplication error of a BRM whose counter starts with an 
arbitrary value. Eq. (2.12) is the error formula resulting when che 
maximum value of the actual output is considered at the points of 
discontinuities. Eq. (2.13) is the companion equation resulting 
when the minimum value of the actual output is considered at these 
points. In this section these error formulas will be analyzed with 
the objective of obtaining error bounds for a BRM with this added 
degree of freedom. We begin by analyzing Eq. (2.12) to obtain the 
maximum positive error of an n stage BRM. 

It is convenient for this discussion to define a vector b 
such that 


b -ll 

b-2 

1 

i 

/y-ii 

y -2 

• 

= M 

I . 

• 

b -nj 


• 

* 

[y-n 


(B.l) 


Theorem B.l For all x^ and xg^ in Eq. (2.12), G < j b.j | ♦ 

i=l 

n 

Moreover, if y..^ = x^ - xg^ then G ^ |b_-jJ where x de- 


i-1 


Sk* 


notes the complement of x 
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The elements of the vector defined by Eq. (B.,l) are 

, 1 (l 1 1 \ 

”-i ■ 2 y-X - \4 y -2 + 8 y -3 + • • • + ^ y -n) 


. 1 fl ^ 1 

b -2 - 2 ^2 ’ (4 y -5 + a y -4 ■*' • 


• + 2 -^h 


, 1 /l .1 

b -k = 2 y -k “ \T y-k-1 + 8 y-k-£ + 


gn-k+l 


r ' n ) 


t>-n ~ 2 y~n 
l/2 > 1/4 + . . . + 
b -k 


It is next noted that each element (x^ - x^) may have only 
three possible values; that is, 0,1, or -1. This may be verified 
by direct computation. 


l/2 n 

then 


< 0 

if 

1 

W 

II 

O 


> 0 

if 

y -k = 1 

(B.2) 


x k 

x Sk 

x k “ x sk 

0 

0 

0 

0 

X 

-1 

1 

0 

1 

1 

1 

0 


(B.3) 


Since | x^ - xgjJ <1 then 


n 


G = (xi - + (x 2 - Xggjb^g + • • • * (x n - *g n )k_ n < ^> (bj 

i=l 


(B.4) 


A sufficient condition for G to attain the upper bound of 
n 

Eq. (B.4), that is, S. IMI. 

ft 


is that: 



13 4 


( *k - * S k ) ‘ - 1 

if 

b_- 1 0 


a +1 

if 

b -k > 0 

(B.5) 


Combining (B.2), (B.3), and (t.o) results in 

y -k [ b -k | x k “ x Sk _ x kj x Skj 

0 <0 ! -1 0 1 

1 i >0 ! +1 1_J 0 ; 

Therefore, y_ k = x k = Xgj, is a sufficient condition for 

»=J> 

i=l 

As a consequence of Theorem B.l the maximum of G; denoted by 

n 

G max , is such that G mx = max 5> ! b _i! • ^ procedure Aich is to be 

7 i -1 n 

followed is to find the value y where the maximum of \ j b_j | 

i=l 

is attained and then evaluating this function. In order to aid this 
analysis the notation b„^ (y) is introduced, where 
y = * • ->y_ n ana b_j(y) denotes the value b_^ for the 

vector (y_- L ,y_ 2 , . . *,y_ n ) . For an n stage BRM, all possible 
b_i(y) values may be obtained by multiplying M with all possible 
values of y. We will call this particular matrix B r . 
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CD 


PQ 



H H H 

I t S 


3 


03 





£ 

03 * • • * 

v-/ 

03 

t 

,0 






S H 

03 


■? 

G 

03 


1 


hIoo hIth 

i t 

nf^ ihIos 

t 


H}^ Hl03 

1 


H 1 03 


Hl03 O • • • 


j* 

pq 



HM< Hi l03 

i 


II 




Three-stage 
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H I CO HN H lc\3 
Hlcvj o 

rolca HN Hloo 
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Hi 03 o o 

tO loo H loo 

I 


riM H f 0*3 O 

S 

HlcO Hl'tf H lev} 

_ 1 I 


H |od Hi-# H icvj 

t i 

H ICVJ O 

! 

Hlcvj O O 


h~ 

!! 

H 

H 

hIoj o o o 

I'jQ 

^ !H HlcO HN* H f 00 

t 

H 

H 

O 

H 

O 

H 

^Ih hN* hIoo o 

H 

O 

O 

i 

O 

H 

H 

10 |H foloo Hh# Hlcvj 

O 

H 

O 

O 

O 

H 

t 1 


^IhhIcvj o o 

I 

*°!h to loo 


m Ih hN< 


•“Mh hIco 

t S 
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Theorem B.2 . For all values of y, |^.^(y)| ° |b^(2 n -y)| 

This result follows hy induction on the number of stages k. 

It was shown as an example for the k = 2 case. Assume it is true 
for the k ** n - 1 case. Then the k » n case is 

/ 1 2 2 n ~ 1 - 1 1 2 11 - 1 - 1 

r — — - . . » s* rr - • « 


®n 


2 n 2 n 


2 n 


2 

0 


2 n 


2 n 


B 


n-x 


®n-l 


where this case is partitioned to show its structure. This theorem 
is obviously true for the first row. The b„^(y) element for the 
n - 1 case is now the b^j^Cy) and the b^^(2 n + y) elements of 
the n case; and the b_i(2 n - y) element of the n - 1 case is 
-ow the b_£_^(2 n “ y) a ^d the b^^ i _^(2 n + 2 n - y) elements of the 
Bjj case. By the induction hypothesis |b_^(y)j = jb„ i (2 n - y) j 
for the n - 1 case. Therefore, these elements fu* the nth cab ; 
yield 

|t>-i-i(y)| - jB-i-i(2 n+1 - y)| >b.7) 


and 


|B - 1 _ 1 ( 2 n - y)| = iKuite 11 + y)l 
Substituting u ® 2 n - y into Eq.. (B.8) then 

l b -i-a (u >l a Ki-l * 2 * 1 * 1 " u )l 


(B.8) 


which completes the proof. 
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Lemma B.l. There exists a y* in the domain 
010... 0 0 < y* <0111. .11 such that 

y_ ki(y*)l > ^ h,j(y)l » 

i" i 


This theorem states that G max is attained in the domain 


010 .. . 00<y<011 ... 11. Asa result of Theorem B.l, 
the search for a point where G mny is attained can be immediately 

restricted to the y domain 0<y < 1 0 0 . . .0 0. Consider B^ 
for these values of y. 


1 1 

4 * * * 2 * 


Bn- 


The vector 



1 

2 


0 % 


0 

0 . . . 

♦ 

» 

0 . . . 


l/2 J results from the y vector 
0 


/ 0 
1 
0 
0 


vector 


1 \ 

0 

0 


0 


can be immediately ruled out. 


'T\ 


The 


The structure Bj, 


Br the above matrix la Vl 
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Because of Theorem B.2, the absolute values of the elements in 
are identical to the absolute values of the elements Bp. More- 
over, since the elements of first row, that is, |b^{y)J, increase 

as y increases then for each column sum to the left of f»l/4i 

1 / 2 | 

0 


there is a column sum to the right which exceeds it Therefore, 
must lie to the right. 

As a result of Theorem B.2 and Lemma B.l, there must be at 
least two values of y where G max is attained. For an n stage 
BRM, we will call the y value corresponding to G^^ on the left 
of y = 1000 . . . 00, and the one on the right of 

y - 100 . . . 00, Rjj. 

Lemma B.2. For an n-stage BRM 

0R n-l < < 0 1 1 . . .11 

This result follows from the proof of Lemma B.l. Since the 
first row jb.^Cy) | increasing then G^^ must lie between the 
right maximum of B^i and the rightmost value of Iiemma B.l. 

Lemma B.5. For an n-stage BRM 


0100. . ♦0<L n <01 1^.2 

Consider Bq for the values of y of Lemma B.l, that is . 
0100. . . 0 < y < 0 1 1 . . .1. 
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1 2 n ~ 2 + 1 2 n " 2 4 - 2 

4 ~ 2 n 2 n 


1 2 n ~ 2 - 1 2 P ~ 2 - 2 

2 gn-l on-1 


Bn-2 



The sums of the absolute value of the first two rows are 


3 Z 11 " 1 + 2 n *" 2 - 1 a 11 " 1 + 2 n ** 2 - 2 

i f on ; on 


• • • 


T’leivf ; re, this is a decreasing sequence. Since by Theorem B.2 and 
Lemma 3.1. B^g attains at a maximum for. at .least' two values, then 
Gma>; must lie between the leftmost value of Lemma B.l and OlL^g.. 
Theorem B.3 . <?cr ah n-stage BPM 

. 0R n-l < ^n < 01 %-2 

..Juls theorem is the combi —on of Lemma B.2 and Lemma 
I h e oren' 1.4. R n and L n are unique and 0R n _y = L^ = l^g • 

This theorem follows immediately from the proofs of 1 -iuma B.2 
and Lemma B.3 by using the principle of strong induct as the 
method of proof; that is, assume it is true for k < n and prove 
for the case k = n. 

The values of y where G^. is attained can be found by 
using Theorem B.4. These values are listed in Table B.l.' 
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TABLE B. 1 VALUES OP y WHERE IS ATTAINED 




<*n-l 

hi 

R n 


BBS 

PWP1 

BWHI 

HBW 




i 

a 1 

2 

II 


01 


3 


Oil 

Oil 


4 


0101 

0101 


5 


01011 

01011 

10101 

6 

010101 

010101 

010101 

101011 


The BRM counter value and starting value corresponding to the 
y values listed in Table B. 1 can be obtained by Theorem B. 1. These 
values are tabulated in Table 2.6. 

An equation for Gyn nY as a function of n may be obtained by 
a procedure similar to that used to obtain E max . In particular, 
using the pattern estaolished above a difference equation may be 
written for n even, and a difference equation may be written for 
n odd. Combining the solutions of these difference - equations, '■ 
may be obtained as a function of n. 


<W n > 


— “ *7 


(B.9) 


9 3 9* 2 n 

Noting the similarity between the rightmost term of Eq. (2. 13)- 
to that of Eq. (2.12), one may establish immediately a minimum error' 
bound when the BfM counter starts with an arbitrary value. /o - 




9 ' 3 g. 2 & 

A tight error bound may be obtained as follows; Expanding 
Eq. (2.13), then equation for ~H may be expressed as 


(B.10) 


H = -xgjjy^ - - C sl )b_ 1 + (Co - C S2 )b. 2 


+ . 


' ( c n “ %h^-nj 

. (B.11) V - 
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x gR identifies the rightmost 1 in the initial value of the 
BRM counter. It may be simply argued that if a binary number has a 
rightmost 1 n position R then its 2 f s complement also has a 1 in 
that position and, moreover, has oeros at all positions j where 
j < R. Therefore, Xg R = Cg R and Cgj = 0 for all j < R. 

Eq. (B.1JL) may be expressed as 

"H = c SRy-R + [^l b -l + C 2 b -2 + • • * + ^R-l b -R+lj 

+ (C R - CgR)b- R + [( C R+1 “ c SR+l) b -F.-l + • • • + (C n -C Sn )b_ n j 

(B.12) 


or alternately as 

-H = Cjb..n + C R b.2 + . . ♦ + C R b_ R + (Cr+i - CsR+l) b -R-l 


+ • • • + (c n - c Sn )b_ n + c. 


SR I 


1 a 1 

2 y-R + 4 y-R-i + * 


*1 

j 

(B.13) 


Therefore 

< | 4 jC2 b -2| + • • • + j ^R b -R | 


+ K°R+1 " c SR+l) b -R-ll + * * • + K c n “ c Sn) b -nl 


+ C, 


SR 


fl 


2 y-R 


i 

4 


+ * y — R— i "** 


(B.14) 


The upper bound of -H is attained when equality is attained 

in Eq. (B.14) . It can be demonstrated by an argument similar to 

that used in Theorem B.l, that the conditions for equality are 

y . = C. for all j 
-0 J 

and 


% = C i 


for all j > R 
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Therefore, Eq. (B.14) may be written as 

-H = (C^C^Cj, . . .,C n )M/C2t 

C 3 


'til 


~ (0,0,0,. . . 0,Cjj + 2,0-j + 2 


> • • • 


Cl 1 

°2 


C3 1 


+ CgR 


\ c r 4 r c R+i + • • * 


(B.15) 


|2 4 

The last two terms in Eq. (B.15) is always nonnegative. 
Moreover, the sum of these two terms is nondecreasing as h de- 
creases. Since by Theorem A.l, = 1 is a condition for maximiz- 
ing the first term of Eq. (B.15), it must also be a condition for 
maximizing Eq. (B.15) itself. 

Eq. (B.ll) may be rewritten with this condition as f ol.l owe 


-H - y -]L + (Cg - 0g 2 ,C 3 - 0 g3 , . . .,C E - Og n )M 


y-2 




(B.16) 


Therefore, the equation for IL^ can be immediately 
written as: 

W n > “ - 1 - Wa- 1 ) “ -X - ^ 


(E.17) 
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Basel on this analysis Cg^, x si^ x ±> an< ^ y_j are equal 
to 1 for The remaining stages are determined such as to 

maximize G for an n - 1 stage BRM. The values of y, C, and 
Cg where H^r, is attained are listed in Table B.3. 


TABLE B.3 C,; y, AM) Cg WHERE H min IS ATTAINED 


; n 

y-i— y-e 


C 6*‘* C 1 

C S6*** C S1 

i 

2 

11 

n 

01 

1 

o 

101 

101 

on 

4 

1011 

1101 

oon 

5 

10101 

10101 

01011 

6 

101011 

110101 

001011 

i 

y«i- * -y-e 

C 6*** C 1 

C S6* * * C S1 

i 

! 2 

ii 

11 

01 

3 

in 

in 

001 

4 

1101 

ion 

0101 

• 5 

non 

non 

00101 

6 

noioi 

loion 

010101 


The BRM counter value and initial value for these H^n 
■values are the 2’s complement of the above numbers. These values 


are tabulated in Table 2.7 



